Listas de estaciones según entidad y estado

Paquetes

library(sf)
library(leaflet)
library(leaflet.extras)
library(mapview)
library(readODS)
library(readxl)
library(parzer)
library(tidyverse)
library(kableExtra)
library(fuzzyjoin)
library(janitor)
source('R/funciones.R')
leaflet_map_view <- . %>% setView(lat = 18.7, lng = -70.3, zoom = 8)

Resultados

En el archivo combinadas_v0.9.xlsx se consolidaron los resultados de los archivos aportados por ONAMET e INDRHI en un único archivo combinado. Este producto derivado es, por lo tanto, una relación de estaciones (en filas) con un conjunto de atributos (columnas) específicos de cada una. Hasta el momento, se han consolidado las de ONAMET (primeras cuatro hojas) e INDRHI (últimas dos hojas).

ONAMET

ONAMET aportó dos archivos, cada uno conteniendo listas de estaciones que entendemos son comunes. El primero contiene una relación de estaciones con indicación de coordenadas (COORDENADAS 2022 DIVISION INSTRUMENTOS.xlsx), mientras que el segundo indica los periodos de medición de cada variable de cada estación (LISTADO PERIODO DE MEDICION DE CADA ESTACION(22)-1.xlsx). El objetivo de este análisis es producir un mapa de las estaciones que se consideran activas, para lo cual necesitamos las coordenadas del primer archivo, y el periodo de medición del segundo.

Para decidir cuáles estaciones se consideran activas, se buscó, mediante expresión regular, la cadena de caracteres “2021” en las columnas que contenían fechas (e.g. columnas PR, TX, etc.) de las listas de estaciones con indicación de periodos de medición. Cada coincidencia señalaba una estación que colecta datos hasta 2022. Se buscó “2021” como equivalente de 2022 porque, como bien explicó ONAMET en el archivo fuente, el año actual (2022) no se hace constar hasta tanto haya transcurrido al completo.

archivo_combinadas <- 'combinadas_v0.91.xlsx'
onamet <- map(1:4, ~ read_xlsx(paste0('fuentes/combinadas/', archivo_combinadas), sheet = .x))
lista_con_periodo <- bind_rows(sapply(2:4, function(x) onamet[[x]], simplify = F)) %>% 
   rename_at(., vars(-'id',-'nombre'), ~ paste('var_', .))
indice_2022 <- sapply(lista_con_periodo %>% select(3:last_col()),
       function(x) grep('2021', x), simplify = F) %>%
  unlist(use.names = F) %>% unique %>% sort
total_2022 <- length(indice_2022)
#Estaciones que colectaban hasta 2022
lista_con_periodo_2022 <- lista_con_periodo[indice_2022,]
#Estaciones que colectaban hasta 2022 con correspondencia en tabla EMC (coordenadas)
lista_con_periodo_2022_coord <- lista_con_periodo_2022 %>% inner_join(onamet[[1]], by='nombre')
anadidas <- which(lista_con_periodo_2022_coord %>% pull(lat) %>% is.na())

En la lista de estaciones meteorológicas convencionales con indicación de coordenadas, se encontraron 83, mientras que la lista de estaciones con indicación de periodo de medición contenía 77 estaciones. De las estaciones con indicación de periodo, se encontraron un total de 36 estaciones registrando datos hasta 2022. De estas, a 32 estaciones se les localizaron sus homólogas en la lista de coordenadas, y a 4 (CABO ENGAÑO, LA VEGA, LOYOLA (SCR ), SANTIAGO) no se les pudo localizar sus correspondientes homólogas con coordenadas. En los casos de LA VEGA y SANTIAGO, la incertidumbre se debió a que hay varias estaciones en dichas ciudades, lo cual hacía difícil determinar con cuáles estaciones hacerlas corresponder. Por lo tanto, en los archivos facilitados, se encontraron 87, de las que 36 se encontraban en estado “activa” y las restantes 51 en estado “inactiva o no reportada”.

Como se puede notar, el mayor desafío que presentaron los datos de ONAMET, fue determinar la correspondencia de las estaciones entre listas debido a la ausencia de un código (clave) común. Si bien cada lista de estaciones asignaba códigos o nombres únicos a las estación, estos no eran consistentes entre listas. Aunque muchas estaciones pudieron emparejarse, en algunos casos hicimos la correspondencia con incertidumbre o, directamente, no se pudo hacer la relación como ya se refirió (ver también comentarios en globos en el archivo combinadas_v0.91.xlsx). Entendemos que si cada estación contase con un identificador único, la correspondencia se hubiese realizado de forma más eficiente y certera.

Generamos una columna de estado, asignando estado “activa” a las que recogen datos hasta 2022. Al mismo tiempo, para elaborar el mapa en secciones posteriores, asignamos coordenadas ad hoc a las estaciones que carecían de ellas. Incluimos la tabla resumen a continuación:

lista_con_periodo_2022_coord_todas <- lista_con_periodo_2022 %>%
  right_join(onamet[[1]], by='nombre') %>% rowwise() %>% 
  mutate(inactiva = all(is.na(c_across(starts_with('var_'))))) %>% 
  mutate(estado = ifelse(
    inactiva,
    'inactiva o no reportada',
    'activa')) %>%
  select(id, nombre_onamet, starts_with('var_'), lat, lon, h_m, estado) %>% 
  mutate(
    lat = case_when(
      nombre_onamet == 'LOYOLA (SCR )' ~ 18.411929,
      nombre_onamet == 'SANTIAGO'  ~ 19.50054,
      nombre_onamet == 'CABO ENGAÑO'  ~ 18.616784,
      nombre_onamet == 'LA VEGA'  ~ 19.22275,
      TRUE ~ lat),
    lon = case_when(
      nombre_onamet == 'LOYOLA (SCR )' ~ -70.112727,
      nombre_onamet == 'SANTIAGO'  ~ -70.69806,
      nombre_onamet == 'CABO ENGAÑO'  ~ -68.325500,
      nombre_onamet == 'LA VEGA'  ~ -70.53108,
      TRUE ~ lon),
    estado = case_when(
      nombre_onamet == 'BANÍ' ~ 'inactiva o no reportada', #Lista de estaciones con coords.
      # nombre_onamet == 'LA VEGA - IATESA' ~ 'inactiva o no reportada', #Lista de estaciones con coords.
      TRUE ~ estado)
    ) %>% 
  select(id, nombre_onamet, lat, lon, h_m, estado) %>% 
  inner_join(lista_con_periodo %>%
              select(nombre, starts_with('var_')) %>% 
              right_join(onamet[[1]] %>% select(nombre, nombre_onamet))
            )
# lista_con_periodo_2022_coord_todas %>% write_csv('estaciones_de_ONAMET_con_indicacion_de_estado.csv')

lista_con_periodo_2022_coord_todas %>%
  arrange(nombre_onamet) %>% 
  kable(booktabs=T) %>%
  kable_styling(latex_options = c("HOLD_position", "scale_down")) %>%
  gsub(' NA ', '', .)

id	nombre_onamet	lat	lon	h_m	estado	nombre	var_ pr	var_ tx	var_ tm	var_ pe	var_ hs	var_ dv	var_ vv	var_ vvmx	var_ nb	var_ hr	var_ ro	var_ tv	var_ ev
	ALTAMIRA	19.67516	-70.83400	272	inactiva o no reportada	ALTAMIRA	1950-2013	1951-82	1951-82
	ANGELINA	19.13138	-70.21977	54	inactiva o no reportada
78466	ARPT. ARROYO BARRIL	19.20050	-69.43144	37	activa	A. BARRIL	1976-88. 91-2021	1976-2021	1976-2021	1980-2021		1976-2021	1995-2021	1976-2021	1980-2021	1977-2021	1976-2021	1976-2006
21199	ARPT. DE LA ROMANA	18.44853	-68.90925	62	activa	LA ROMANA	1919, 1923-27, 1931-80, 2006-2021	1924-27, 1931-80, 2006-2021	1924-27, 1931-80, 2006-2021	2006-2021	2006-14	2006-2021	2006-2021	2006-2021	2006-2021	2006-2021	2006-2021	2006-2021
78485	ARPT. DE LAS AMÉRICAS	18.43306	-69.67959	7	activa	LAS AMERICAS	1957-2021	1960-2021	1960-2021	1957-2021		1957-2021	1957-2021	1956-13	1956-2021	1961-2021	1968-2021	1968-2021
	ARPT. DE PUNTA CANA	18.56400	-68.35941	7	inactiva o no reportada
	ARPT. DEL CIBAO	19.40305	-70.59781	170	inactiva o no reportada
	ARPT. EL CATEY	19.26700	-69.73366	4	inactiva o no reportada
78457	ARPT. GREGORIO LUPERÓN	19.75416	-70.56316	16	activa	A. LA UNION	1977-2021	1977-2021	1977-2021	1976-2021	1977-83	1977-2021	1977-2021	1977-2021	1977-2021	1977-2021	1977-2014	1977-2021
	ARPT. LA ISABELA	18.57696	-69.98158	22	inactiva o no reportada
78482	ARPT. MARÍA MONTEZ	18.24861	-71.12288	6	activa	BARAHONA	1931-2021	1951-95-2021	1951-2021	1975-2021	1975-88	1951-2021	1961-2021	1960-2021	1957-2021	1960-2021	1968-2021	1968-2021	1976-95
	ARROYO LORO	18.81403	-71.27818	419	inactiva o no reportada
	BANÍ	18.27867	-70.31014	57	inactiva o no reportada	BANI	1936-2008	1938-93	1938-93
	BÁNICA	19.08156	-71.69836	285	inactiva o no reportada
	BASE AÉREA DE SAN ISIDRO	18.49944	-69.76989	33	inactiva o no reportada
78473	BAYAGUANA	18.74216	-69.63083	53	activa	BAYAGUANA	1938-2021	1958-2021	1958-2021	1976-2021		1952-2021		1987-2021	1977-2021	1976-2021	1975-2021	1975-2021
	BOHECHÍO	18.77747	-70.99150	478	inactiva o no reportada
	BONAO	18.93317	-70.40500	175	inactiva o no reportada	BONAO	1939-96	1937-97	1937-97		1975-87
78478	CABO ENGAÑO	18.61678	-68.32550		activa	CABO ENGAÑO	1951-2021	1951-2021	1951-2021	1976-2021	1976-77	1951-2021	1951-2021	1976-2021	1951-2021	1952-2021	1968-2021	1968-2021
	CABRAL	18.24784	-71.21920	26	inactiva o no reportada	CABRAL	1939-2003	1958-81	1958-81		1975-88
78464	CABRERA	19.64436	-69.90633	18	activa	CABRERA	1938-2003,06-07,09-2021	1960-83,06-07,096-2021	1960-83,06-07,096-2021			1981 -94,13-2021	1981 -94,13-2021	1975-90		1976-1999	1975-1986	1975-86
22584	CONSTANZA	18.91173	-70.71619	1202	activa	CONSTANZA	1931-91, 94-2013,2021	1931-2021	1931-2021				1966-91			1966-2007
	COTUÍ	19.04825	-70.15273	79	inactiva o no reportada	COTUI	1938-2003	1938-40, 58-87	1938-87
	DAJABÓN	19.54325	-71.70436	44	inactiva o no reportada	DAJABON	1931-2003	1959-2003	1959-2003				1959-68
	DUVERGÉ	18.37666	-71.51814	1	inactiva o no reportada	DUVERGE	39-50, 60-85, 94-2004	1961-85	1961-85
	EL CERCADO	18.72556	-71.51884	739	inactiva o no reportada	EL CERCADO	1939-2003	1958-96	1958-96
	EL SEIBO	18.77438	-69.03995	94	inactiva o no reportada	EL SEYBO	1945-92	1945-92	1945-92
	ELÍAS PIÑA	18.87600	-71.70133	402	inactiva o no reportada	ELIAS PIÑA	1938-97	1952-87	1952-87
	ENRIQUILLO	17.89165	-71.24173	14	inactiva o no reportada	ENRIQUILLO	1939-92	1967-86	1967-86
	GASPAR HERNÁNDEZ	19.62383	-70.27311	13	inactiva o no reportada	GASPAR HDEZ.	1938-2013	1964-87	1964-87
	GURABO	19.49794	-70.66713	244	inactiva o no reportada
	HATO MAYOR	18.77005	-69.25557	104	inactiva o no reportada	HATO MAYOR	1934-90	1951-89	1951-89		1983-95
	HIGÜEY	18.62045	-68.70098	86	inactiva o no reportada	HIGUEY	1934-87, 95-96	1936-87	1936-87
22754	HONDO VALLE	18.71213	-71.69364	882	activa	HONDO VALLE	1953-2021	1953-1998,13-2021	1953-1998,13-2021
	IMBERT	19.75626	-70.83042	128	inactiva o no reportada	IMBERT	1939-2003	1951-77	1951-77
23523	JARABACOA	19.12551	-70.60572	570	activa	JARABACOA	1931-2021	1932-82,86-01,11-2021	1932-82,86-01,11-2021		1976-95
	JARDÍN BOTÁNICO NACIONAL	18.49489	-69.95293	48	inactiva o no reportada
78480	JIMANÍ	18.49278	-71.85301	45	activa	JIMANI	1948-2021	1961-2021	1961-2021	1985-2021		1948-2021	1979-2021	1976-2021	1976-2021	1976-2013	1977-2013	1976-2021
	JUMA BONAO	18.90047	-70.38549	181	inactiva o no reportada
23551	LA VEGA	19.22275	-70.53108		activa	LA VEGA	1924-2006,09-2021	1941-92, 94-06,09-2021	1941-92, 94-06,09-2021			1941-06			1941-99,11-2021	1977-99,11-2021
	LA VEGA - IATESA	19.25230	-70.55246	133	inactiva o no reportada
	LA VEGA - M. A.	19.19920	-70.49850	82	inactiva o no reportada
22346	LA VICTORIA	18.59779	-69.83956	17	activa	LA VICTORIA	1938-89,10-13,2021	1958-89,95, 97,10-13,2021	1958-89,95, 97,10-13,2021
	LAS MATAS DE FARFÁN	18.87478	-71.52819	430	inactiva o no reportada	LAS M. DE FARFAN	1948-2003	1958-2000	1958-2000
	LAS MATAS DE SANTA CRUZ	19.65356	-71.50686	44	inactiva o no reportada
	LIMÓN DEL YUNA	19.11985	-69.78488	11	inactiva o no reportada
	LOMA DE CABRERA	19.41925	-71.61228	227	inactiva o no reportada
22239	LOS LLANOS	18.62499	-69.49096	41	activa	LOS LLANOS	1940-03,2013-2021	1977-86, 95, 97	1977-86, 95, 97
	LOYOLA (SCR )	18.41193	-70.11273		activa	LOYOLA (SCR )	1981-2021	1981-2021	1981-2021						1981-99	1981-99			1981-99
	LUPERÓN	19.88717	-70.96424	10	inactiva o no reportada	LUPERON	1949-2001	1949-51, 93-98	1949-51, 93-98
	MAO	19.58833	-71.04966	65	inactiva o no reportada	MAO	1938-2007	1947-97	1947-97
	MICHES	18.98099	-69.04221	14	inactiva o no reportada	MICHES	1931-81, 96-97	1957-82	1957-82
23591	MOCA	19.38811	-70.53213	179	activa	MOCA	1931-2021	38-87, 95-2021	1938-87,95-2021
23693	MONCIÓN	19.40988	-71.16163	382	activa	MONCION	1931-2007,2014-2021	31-42, 51-96,14-2021	1931-42, 51-96,14-2021
78451	MONTE CRISTI	19.84988	-71.65350	8	activa	MONTE CRISTY	1933-2021	1951-2021	1951-2021	1976-95,2007-2021	1976-88		1976-2021		1976-2021	1976-2021	1976-2021	1976-2021
	MONTE PLATA	18.80243	-69.78160	47	inactiva o no reportada	MONTE PLATA	1938-91	1952-87, 95, 87	1952-87, 95, 87
	NAVARRETE	19.56062	-70.86739	137	inactiva o no reportada
	NEIBA	18.48347	-71.41889	32	inactiva o no reportada	NEYBA	1943-88, 95-03	1950-92	1950-92
20668	OVIEDO	17.80205	-71.40403	35	activa	OVIEDO	1963-87, 94-03,2009-2021	1965-85,13-2021	1965-85,13,2021
22559	PADRE LAS CASAS	18.73026	-70.94115	505	activa	P. LAS CASAS	1937-2014	1950-82,13-2021	1950-82,13-2021
	PARQUE MIRADOR NORTE	18.52763	-69.93112	22	inactiva o no reportada
	PARQUE MIRADOR SUR	18.45298	-69.93416	44	inactiva o no reportada
	PERALTA	18.57806	-70.75596	556	activa	PERALTA	1939-2021
21626	POLO	18.08746	-71.27231	713	activa	POLO	1948-2021	1951-2021	1951-2021		1978-84				1976-99	1966-99
	RADIOSONDA	18.47416	-69.87016	16	inactiva o no reportada
22449	RANCHO ARRIBA	18.69438	-70.44973	681	activa	RANCHO ARRIBA	1939-2021	1961-81,2002,2021	1961-81,2002,2021			1971-00			1976-00
	RESTAURACIÓN	19.31497	-71.68913	579	inactiva o no reportada	RESTAURACION	1939-98	1958-87	1958-87
24421	RÍO SAN JUAN	19.63251	-70.07475	17	activa	RIO SAN JUAN	1948-2021	1961-83,14-2021	1961-83,14-2021
78467	SABANA DE LA MAR	19.05267	-69.38883	11	activa	S. DE LA MAR	1939-2021	1940-2021	1940-2021	1976-2021	1976-88	1951-2021	1957-2021	1951-2021	1951-2021	1957-2021	1968-2014	1969-2021
	SABANA GRANDE DE BOYÁ	18.94384	-69.80099	280	inactiva o no reportada
23489	SALCEDO	19.37531	-70.41498	193	activa	SALCEDO	1931-2021	1939-2011	1939-2011
23256	SAMANÁ	19.20438	-69.33286	9	activa	SAMANA	1931-2021	1951-87,10-2021	1951-87,10-2021
	SAN CRISTÓBAL	18.40710	-70.09011	23	inactiva o no reportada	SAN CRISTOBAL	1934-2005	1942-94	1942-94		1976-88	1949-84	41-75, 79, 88	1941-79	1951-2006	1942-2006			1976-87
22510	SAN JOSÉ DE OCOA	18.54200	-70.50950	469	activa	SAN J. OCOA	1931-2002,10-2021	1931-86,10-2021	1931-86,10-2021
21193	SAN RAFAEL DEL YUMA	18.42610	-68.67416	51	activa	SAN R. YUMA	1943-2021	1952-2021	1952-2021
23352	SÁNCHEZ	19.23314	-69.61311	32	activa	SANCHEZ	1931-2021	1939-84,12-2021	1939-84,12-2021
78460	SANTIAGO	19.50054	-70.69806		activa	SANTIAGO	1931-2021	1941-2021	1941-2021	1976-2021	1975-86	1943-2021	1959-2021	1959-2021	1941-2021	1959-2021	1968-2021	1968-2021	1976-95
	SANTIAGO RODRÍGUEZ	19.47712	-71.33899	130	inactiva o no reportada	SANTIAGO RGUEZ.	38-40, 43,47-12	39-43, 59-65,77-82	1939-82
78486	SANTO DOMINGO ESTE	18.47339	-69.87053	14	activa	SANTO DOMINGO	1931-2021	1936-2021	1936-2021	1975-2021	1984-95	1947-2021	1951-14	1951-2021	1931-2021	1951-2021	1968-2021	1968-2021	1960-2021
	TÁBARA ABAJO	18.46438	-70.86903	126	inactiva o no reportada
	TAMBORIL	19.48342	-70.61141	257	inactiva o no reportada
22443	VILLA ALTAGRACIA	18.67582	-70.17060	183	activa	VILLA ALTAG.	1938-2021	1951-95,13-2021	1951-95,13-2021
	VILLA ISABELA	19.82111	-71.05694	32	inactiva o no reportada
	VILLA LOS ALMÁCIGOS	19.41080	-71.44328	214	inactiva o no reportada
23358	VILLA RIVA	19.17993	-69.91483	25	activa	VILLA RIVA	39-40, 52-80,87-2021	39-40, 52-73, 77-81,14-2021	39-40, 52-73, 77-81,14-2021
24659	VILLA VÁSQUEZ	19.74207	-71.45730	32	activa	VILLA VASQUEZ	1939-2021	1951-1986,13-2021	1951-1986,13-2021
	YAMASÁ	18.77271	-70.02175	84	inactiva o no reportada	YAMASA	1938-2002	1952-83, 95, 97	1952-83, 95, 97

La siguiente tabla muestra un resumen por estado:

onamet_estado <- lista_con_periodo_2022_coord_todas %>%
  group_by(estado) %>%
  summarise(n = n())
onamet_estado %>% 
  kable(booktabs=T) %>%
  kable_styling(
    full_width = F,
    latex_options = c("HOLD_position", "scale_down")) %>%
  gsub(' NA ', '', .)

estado	n
activa	36
inactiva o no reportada	51

El mapa de las estaciones de ONAMET, según el estado actual de las mismas, se muestra a continuación. Se recuerda que el mapa representa, cartográficamente, lo contenido en la lista de estaciones y los correspondientes periodos reportados de operación por ONAMET.

lista_con_periodo_2022_coord_todas_sf <- st_as_sf(lista_con_periodo_2022_coord_todas,
                                                  coords = c('lon', 'lat'), crs = 4326)

fpal_estado <- colorFactor(
  palette = RColorBrewer::brewer.pal(length(unique(lista_con_periodo_2022_coord_todas_sf$estado)), 'Set1'),
  domain = unique(lista_con_periodo_2022_coord_todas_sf$estado), reverse = T)
leaflet(lista_con_periodo_2022_coord_todas_sf) %>%
  addCircleMarkers(
    radius = 5, label = ~ nombre_onamet, group = ~ estado,
    color = ~ fpal_estado(estado),
    stroke = F, fillOpacity = 1
  ) %>%
  addLegend(pal = fpal_estado, values = ~ estado, opacity = 1,
            title = "Estaciones<br>Meteorológicas<br>Convencionales<br>ONAMET<br>Estado") %>% 
  addTiles(group = 'OSM') %>%
  addProviderTiles("Esri.NatGeoWorldMap", group="ESRI Mapa") %>%
  addProviderTiles("Esri.WorldImagery", group="ESRI Imagen") %>%
  addProviderTiles("CartoDB.Positron", group= "CartoDB") %>%
  addLayersControl(
    baseGroups = c("ESRI Imagen", "CartoDB", "OSM", "ESRI Mapa"),    
    overlayGroups = ~ estado, position = 'bottomright',
    options = layersControlOptions(collapsed = FALSE)) %>% 
  leaflet_map_view %>% 
  addFullscreenControl()

# Exportar
ruta_gpkg_est_climaticas_onamet <- 'out/con_indicacion_estatus_onamet.gpkg'
lista_con_periodo_2022_coord_todas_sf %>%
  st_write(ruta_gpkg_est_climaticas_onamet, delete_dsn = T, quiet = T, verbose = F)
ruta_csv_est_climaticas_onamet <- 'out/con_indicacion_estatus_onamet.csv'
lista_con_periodo_2022_coord_todas %>%
  select(ID = id, Nombre = nombre_onamet, Estado = estado, Longitud = lon, Latitud = lat) %>%
  write_csv(ruta_csv_est_climaticas_onamet)

NOTA. Se han generado archivos GeoPackage (out/con_indicacion_estatus_onamet.gpkg) y CSV (out/con_indicacion_estatus_onamet.csv), de la lista de estaciones climáticas de ONAMET, depuradas según el procedimiento explicado en el texto.

INDRHI

En el caso de INDRHI se utilizó un único archivo como fuente (Informe Final inventario estaciones hidrometeorológicas INDRHI, Rep. Dom..pdf), el cual fue elaborado por personal técnico de la institución en 2019, con el objetivo de documentar el estado de sus estaciones hidroclimáticas. El documento nos fue entregado en formato PDF, por lo que hubo que extraer la información programáticamente. Las informaciones sobre el estado de las estaciones fueron complementadas con comunicaciones directas con personal del INDRHI.

indrhi <- map(5:6, ~ read_xlsx(paste0('fuentes/combinadas/', archivo_combinadas), sheet = .x))

Las estaciones reportadas en el informe son de dos tipos generales: hidrométricas y climáticas. De las primeras, las hidrométricas, hay un total de 175 estaciones, mientras que las climáticas representan un total de 61 estaciones. En ambos casos, INDRHI clasificó las estaciones según “estado” en tres categorías: Bueno, Regular, Malo.

Estaciones hidrométricas

indrhi_hidrometricas_resumen_estado <- indrhi[[1]] %>% group_by(Estado) %>%
  mutate(Estado = factor(Estado, levels=c('Bueno', 'Regular', 'Malo'))) %>%
  count()
indrhi_hidrometricas_resumen_estado %>%
  kable(booktabs=T) %>%
  kable_styling(
    full_width = F,
    latex_options = c("HOLD_position", "scale_down")) %>%
  gsub(' NA ', '', .)

Estado	n
Bueno	30
Regular	24
Malo	121

Al momento de redactarse el presente documento, se nos informó que el estado respecto de la línea base de 2019 se ha mantenido bastante similar. En el referido informe, de un total de 175 estaciones hidrométricas, 121 se encontraban en la categoría “Malo”, 24 en estado “Regular” y 30 en la categoría “Bueno”. Las estaciones de la categoría “Malo” se encontraban fuera operación al momento de realizarse el referido informe, y presentaban daños significativos, por lo que su restablecimiento podría resultar complicado. En el caso de las estaciones en estado “Regular”, aunque no se encontraban en operación al momento de realizarse el levantamiento en 2019, su restablecimiento requería una inversión muy pequeña. Finalmente, las estaciones en estado “Bueno” se encontraban en operación normal.

A modo de actualización, según Israel Acosta, encargado del departamento de hidrología de INDRHI, se nos confirmó que el número de estaciones según cada categoría de estado, ha cambiado muy poco en años recientes. Añadió que se reportaron algunas incidencias en años recientes, relacionadas con la estabilidad del personal responsable en terreno. Acosta destacó además que una buena parte de las 24 estaciones en estado regular, continúan siendo recuperables actualmente con una mínima inversión.

En cuanto a la localización de las estaciones, la tabla de las hidrométricas aportada por el INDRHI cuenta con dos pares columnas de coordenadas bajo las denominaciones Latitud, Longitud, Latitud Corregida, Longitud Corregida. El primer par contiene valores de lo que parecen ser coordenadas de falso norte y falso este en la proyección UTM (no se especifica datum, asumiremos WGS84), por lo que no se trata de coordenadas latitud/longitud. El segundo par, con el mismo nombre y sufijo “Corregida”, sugiere una mejora en la precisión, y sus coordenadas (mayoritariamente) están expresadas en sistema sexagesimal de latitud y longitud, aunque también hay algunas UTM entremezcladas.

Para garantizar la consistencia de las coordenadas, escribimos una función, basada en expresiones regulares, que eligía de preferencia el par de coordenadas del grupo “Corregida”; si este par no estaba disponible, entonces elegía el par “normal”. Destacamos también que algunas estaciones no disponían de coordenadas de ninguna clase (e.g. estaciones N° 1, 105, 106, 116), y otras registraban valores claramente erróneos y no atribuibles a ningún sistema de referencia (e.g. estaciones N° 118, 119, 120). Nuestra función excluyó dichas estaciones del mapa y de los análisis posteriores. Todas las coordenadas fueron transformadas a sistema de coordenadas geográficas, WGS84 (EPSG:4326). La tabla exhaustiva de las estaciones, conteniendo las coordenadas consolidadas y depuradas, se muestra a continuación.

coord_limpias <- bind_cols(lapply(
  c('Latitud', 'Longitud', 'Latitud Corregida', 'Longitud Corregida'),
  function(x)
    limpiar_coord(
      mi_vector = indrhi[[1]][,x, drop=T],
      sufijo = ifelse(grepl('Corregida', x), 'correg', 'normal'))))
coord_limpias$utm_x_consolidadas <- with(
  data = coord_limpias,
  expr = ifelse(is.na(utm_x_correg), utm_x_normal, utm_x_correg))
coord_limpias$utm_y_consolidadas <- with(
  data = coord_limpias,
  expr = ifelse(is.na(utm_y_correg), utm_y_normal, utm_y_correg))
indrhi_hidrometricas_coord_limpias <- cbind(indrhi[[1]], coord_limpias)

indrhi_hidrometricas_final <- indrhi_hidrometricas_coord_limpias %>% 
  select(`N°`, utm_x_consolidadas, utm_y_consolidadas) %>% 
  na.omit %>% 
  st_as_sf(coords = c('utm_x_consolidadas', 'utm_y_consolidadas'), crs = 32619) %>% 
  st_transform(4326) %>% 
  mutate(lon_dd_consolidadas = unlist(map(.$geometry, 1)),
         lat_dd_consolidadas = unlist(map(.$geometry, 2))) %>% 
  st_drop_geometry() %>% 
  right_join(indrhi_hidrometricas_coord_limpias, by = 'N°') %>% 
  mutate(lon_dd_consolidadas = ifelse(is.na(lon_dd_consolidadas), lon_dd_correg, lon_dd_consolidadas),
         lat_dd_consolidadas = ifelse(is.na(lat_dd_consolidadas), lat_dd_correg, lat_dd_consolidadas)) %>% 
  arrange(`N°`) %>% 
  select(-(utm_x_normal:utm_y_consolidadas), -(Latitud:`Longitud Corregida`)) %>% 
  relocate(matches('consolidadas'), .after = `Sub-cuenca`)
indrhi_hidrometricas_final %>% 
  kable(booktabs=T) %>%
  kable_styling(latex_options = c("HOLD_position", "scale_down")) %>%
  gsub(' NA ', '', .)

N°	Código	Nombre de la estación	Cuenca	Sub-cuenca	lon_dd_consolidadas	lat_dd_consolidadas	Altitud (m)	Mira - Periodo de registro	Estado	Transmisión de datos	Existe Curva de Gasto	Observación
1		Isabela	Bajabonico	Puerto plata					Bueno	Correo fisico: planilla en papel, mensual	no	La garita esta en mala condiciones, el Oservador esta invalido
2	311001	Batey Excabacion	Higuamo		-69.40147	18.63061	16	2015-2017	Malo	Correo físico: planilla en papel, mensual	Incomp.	el lugar se encuentra con maleza que fueron cortadas, y le falta la mira #3,4,6,7,8,9
3	532001	Puerto Escondido	Lago Enriquillo	Enriquillo/Río Las Damas	-71.58066	18.32454	425	1955-2019	Bueno	Correo físico: planilla en papel, mensual	si	Falta de mantenimiento
4		El Millo	Lago Enriquillo	Río Majagual/Enriquillo	-71.34922	18.54715	80	1984-2002	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
5	230001	Anamuya	Rio Anamuya	Anamuya	-68.64482	18.69947	68	2015-2017	Malo	Telemetría Satelital	Incomp.	La estación hidrométrica solo están las bases de las miras, no tienen limnimetro.
6	45001	Rinconcito	Rio Artibonito	Artibonito/Río Macasias	-71.76762	18.95801	268	1955-2018	Bueno	Correo fisico: planilla en papel, mensual	si	Instrumentos descalibrados, Irresponsabilidad de los lectores y poca supervision
7		Sonador	Rio Artibonito	Artibonito/Río Sonador	-71.57687	18.70537	760	1972-2013	Bueno	Correo físico: planilla en papel, mensual	si	No tiene lector de mira, fue pensionado: La mira esta en perfecta condición.
8	543001	Puertecito	Rio Artibonito	Artibonito/Río Macasias	-71.50987	18.80093	605	1955-2017	Bueno	Correo físico: planilla en papel, mensual	si	Aun se puede considerar de utilidad, porque aun están funcionando y reportando
9	542001	Cajuilito	Rio Artibonito	Artibonito/Río Tocino	-71.59966	19.08263	456	1978-2019	Bueno	Correo físico: planilla en papel, mensual	si	Aun se puede considerar de utilidad, porque aun están funcionando y reportando
10	543002	Ranchitos	Rio Artibonito	Artibonito/Río Macasias	-71.61729	18.93246	380	1955-2019	Malo	Correo físico: planilla en papel, mensual	si	La Sección de mira hay que cambiarla porque cambio sus cauces. Datos no son reales.
11	543101	Pozo Hondo	Rio Artibonito	Artibonito/Río Yacahueque	-71.48591	19.01005	315	1964-2019	Regular	Correo físico: planilla en papel, mensual	si	El muro esta en mala condiciones, La lectora de mira esta enferama
12	540002	Pedro Santana	Rio Artibonito	Río Artibonito	-71.69024	19.10257	278	1956-2019	Regular	Correo físico: planilla en papel, mensual	si	La mira 3 no tene regla
13		El Estrecho	Rio Bajabonico	Bajabonico	-70.94147	19.80553	40	1956-1965	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
14	60001	Imbert	Rio Bajabonico	Bajabonico	-70.82714	19.29927	123	1955-1993	Regular	Correo físico: planilla en papel, mensual	si	Faltan las miras 0,1 y 2. Además no tiene observador.
15	400001	El Recodo	Rio Bani	Río Baní	-70.33945	18.37361	230	1979-2006	Malo	Telemetría Satelital	si	La estacion ya no existe por crecida del rio bani.
16		Los Jengibres	Rio Boba	Boba/Boba	-70.04242	19.44386	17	1968-2013	Malo	Correo físico: planilla en papel, mensual	si	Las miras se deslizaron.
17		La Pinta	Rio Chacuey	Chacuey/Chacuey	-71.55150	19.64413	35	1964-1978	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
18		Manabao	Rio Chacuey	Yaque del Norte	-70.78760	19.07700	865	1963-2001	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
19		La Espensa	Rio Chacuey	Chacuey/Chacuey	-71.55428	19.57802	80	1977-1993	Regular	Correo fisico: planilla en papel, mensual	si	Le falta la mira 1
20	260001	Santa Lucía	Rio Chavón	Río Chavón	-68.93690	18.68206	77	1956-2001	Malo	Telemetría Satelital	si	Es necesario reubicar la estación, la sección esta en malas condiciones. La caseta no tiene limnigrafo, las miras fueron derribadas por crecida del rio.
21	240001	El Mamey	Rio Haina	Río Duey	-68.75111	18.66878	0	1968-2019	Bueno	Correo fisico: planilla en papel, mensual	si	Aun se puede considerar de utilidad, porque aun están funcionando y reportando
22	340004	Los Corozos	Rio Haina	Río Haina	-70.12344	18.52466	63	1982-2007	Malo	Telemetría Satelital	si	La estacion no existe por crecida del Rio Haina, el observador fallecio.
23	10001	Don Miguel	Rio Masacre	Masacre/Masacre	-71.68351	19.49365	83	1955-2019	Bueno	Correo físico: planilla en papel, mensual	si	Aun se puede considerar de utilidad, porque aun están funcionando y reportando
24		Carbonera	Rio Masacre	Masacre/Guajabo	-71.69400	19.59052	30	1964-1967	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
25		La Aduana	Rio Masacre	Masacre/Masacre	-71.71141	19.54847	42	1980-1999	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
26	160003	Cinta Negra	Rio Nagua	Nagua/Nagua	-69.96343	19.28035	39	1980-2018	Malo	Telemetría Satelital	si	Este observador abandono y dejo un sustituto, pero ese sustituto también se fue hace 2 años. Por la apariencia de la caseta, puede ser que este limnigrafo no este funcionando hace mucho tiempo.
27		Villa Nizao	Rio Nizaito	Nizaito	-71.18576	18.01662	121	1955-1998	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
28	380011	Bocaina	Rio Nizao	Río Nizao	-70.47316	18.67280	662	1980-2019	Malo	Telemetría Satelital	si	Existen dos casetas una con limnigrafo no funciona, y la otra para telemétrica no tiene ningún equipo por vandalismo. Estación hidrométrica con 4 bases. Solo tiene las miras 0 y 1, faltan las 2 y 3.
29		Jiguey	Rio Nizao	Rio Nizao	-70.37745	18.54599	573	2015-2017	Regular	Telemetría Satelital	Incomp.	Caseta no tiene candado ni batería.
30		Aguacate	Rio Nizao	Rio Nizao	-70.32149	18.49588	347	2015-2017	Regular	Telemetría Satelital	Incomp.	Caseta sin candado no tiene batería.
31	332321	El Salto	Rio Ozama	Comatillo	-69.60801	18.83988	144	2015-2017	Bueno	Telemetría Satelital	Incomp.	Sección de miras en buenas condiciones. No tiene observador(Cancelado)
32	333101	El Higuero	Rio Ozama	Río Higuero	-69.99933	18.61658	49	1959-1985	Malo	Telemetría Satelital	si	La estacion hidrometrica no existe por crecida del Rio Higuero, es necesario reinstalar.
33	330001	Don Juan	Rio Ozama	Rio Ozama	-69.95100	18.81812	55	2015-2017	Malo	Telemetría Satelital	Incomp.	El rio cambio de cauce a unos 30 metros de la estación hidrométrica, es necesario reubicar la estación.
34	332101	Los Hidalgos	Rio Ozama	Rio Sabana	-69.66440	18.87713	178	2015-2017	Malo	Telemetría Satelital	Incomp.	Sección de miras en condiciones precaria, se recomienda reubicar la estación. No tiene observador (Falleció)
35	331001	El Cacique	Rio Ozama	Río Savita	-69.86123	18.81733	37	1955-2014	Regular	Correo fisico: planilla en papel, mensual	si	Hay que Reinstalar la mira 1,2. estan Dañadas(observador pencionado )
36		Valle de Samaná	Rio San Juan	San Juan/San Juan	-69.33747	19.26323	14	1980-2013	Malo	Correo físico: planilla en papel, mensual	si	Instrumentos descalibrados, Irresponsabilidad de los lectores y poca supervisión
37		San Rafael I	Rio San Rafael	Río San Rafael	-71.13787	18.02911	0	1956-2013	Malo	Correo físico: planilla en papel, mensual	si	Instrumentos descalibrados, Irresponsabilidad de los lectores y poca supervisión
38		San Rafael II	Rio San Rafael	Río San Rafael	-71.14120	18.03241	185	1980-2013	Malo	Correo físico: planilla en papel, mensual	si	Instrumentos descalibrados, Irresponsabilidad de los lectores y poca supervisión
39		San Rafael	Rio San Rafael	Río San Rafael/Majagual	-71.14116	18.03042	156	1980-2013	Malo	Correo físico: planilla en papel, mensual	si	Instrumentos descalibrados, Irresponsabilidad de los lectores y poca supervisión
40	300002	Paso al Medio	Rio Soco	Río Soco	-69.17193	18.56742	15	1957-2018	Malo	Correo físico: planilla en papel, mensual	si	hay que re instalar la mira 1,2,3,4 se la llevo la crecida del rio.
41	47002	Santa Cruz	Rio Yaque del Norte	Y.del B./Maguaca	-71.52177	19.62386	47	1964-1978	Bueno	Correo fisico: planilla en papel, mensual	si	Falta la mira #5
42	45501	Puente Guayubin	Rio Yaque del Norte	Yaque del Norte	-71.39816	19.66163	25	1960-1966	Bueno	Correo fisico: planilla en papel, mensual	si	No hay observador
43	42006	Sabaneta	Rio Yaque del Norte	Y. del N./Bao	-70.97864	19.54565	767	1980-2018	Bueno	Correo fisico: planilla en papel, mensual	si	Falta la mira #0 , observador con edad avanzada
44	42801	La Fortaleza	Rio Yaque del Norte	Y.del N/Janico	-70.83344	19.04703	407	1982-2019	Bueno	Correo fisico: planilla en papel, mensual	si	Solo le falta la #0
45	40008	Puente San Rafael	Rio Yaque del Norte	Yaque del Norte	-71.05942	19.58971	62	1958-2017	Bueno	Correo fisico: planilla en papel, mensual	si	Mira 5 esta Gastada
46	40020	Peña Ranchadero	Rio Yaque del Norte	Yaque del Norte	-71.32977	19.65942	49	1977-2018	Bueno	Correo fisico: planilla en papel, mensual	si	Aun se puede considerar de utilidad, porque aun están funcionando y reportando
47	43001	Inoa	Rio Yaque del Norte	Y. del N/Inoa	-70.95897	19.32357	331	1967-1995	Bueno	Correo físico: planilla en papel, mensual	si	Quien hace la observación es la sobrina.
48	40017	Boma	Rio Yaque del Norte	Yaque del Norte	-70.66224	19.16861	450	1971-2017	Bueno	Correo físico: planilla en papel, mensual	si	Aun se puede considerar de utilidad, porque aun están funcionando y reportando
49	40007	Jinamagao	Rio Yaque del Norte	Yaque del Norte	-70.98202	19.54277		1960-2019	Bueno	Correo físico: planilla en papel, mensual	si	Aun se puede considerar de utilidad, porque aun están funcionando y reportando
50	41202	Jarabacoa	Rio Yaque del Norte	Y. del N./Arroyo Cercado	-70.67840	19.08357	510	1964-1967	Malo	Correo fisico: planilla en papel, mensual	si	Esta estacion hay que reacerla completa
51	40010	Palo verde	Rio Yaque del Norte	Yaque del Norte	-71.55546	19.75778	11	1959-1995	Malo	Correo fisico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estacion, debe ser recontruida. Tiene un observador nombrado
52	41740	Mata Grande	Rio Yaque del Norte	Y. del N./Antonsape Bueno	-70.99611	19.19532	910	1981-2019	Malo	Correo físico: planilla en papel, mensual	si	Esta estación se encuentra en mal estado.
53	45301	Cana Chapetón	Rio Yaque del Norte	Y. Del N./Cana	-71.25967	19.59789	65	1959-1966	Malo	Correo físico: planilla en papel, mensual	si	El rio arrastro la mira.
54	42203	Bejucal	Rio Yaque del Norte	del N./Guajuma	-70.73507	19.18246	707	1981-2000	Malo	Correo físico: planilla en papel, mensual	si	Estacion destruida totalmente.
55		Rincón	Rio Yaque del Norte	Y. Del N./Guayubín	-71.38663	19.52613	128	1964-1995	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
56		La Antona	Rio Yaque del Norte	Y. Del N./Guayubín	-71.40289	19.63025	38	1955-1966	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
57		Inaje	Rio Yaque del Norte	Y. Del N./Inaje	-71.46927	19.42135	162	1960-1961	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
58		Jicomé	Rio Yaque del Norte	Y. Del N./Jicomé	-70.94842	19.62802	121	1957-1958	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
59		La Gorra	Rio Yaque del Norte	Y. Del N./Maguaca	-71.49482	19.51830	112	1977-1977	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
60		Bulla	Rio Yaque del Norte	Y. Del N./Mao	-71.07707	19.41982	295	1967-2001	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
61		Chorrera	Rio Yaque del Norte	Y. Del N./Mao	-71.08481	19.46413	108	1957-1967	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
62		Martínez	Rio Yaque del Norte	Y. Del N./Mao	-71.07953	19.47052	103	1956-1965	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
63		Mata de Jobo	Rio Yaque del Norte	Y. Del N./Yaguajai	-71.33121	19.46497	108	1955-1965	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
64		Manabao	Rio Yaque del Norte	Y. del N./Arroyo Grande	-70.78987	19.06973	901	1955-1979	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
65		Paso Bajito	Rio Yaque del Norte	Y. Del N./Baiguate	-70.66399	19.01495	982	1981-2013	Malo	Correo físico: planilla en papel, mensual	si	Instrumentos descalibrados, Irresponsabilidad de los lectores y poca supervisión
66		Bao	Rio Yaque del Norte	Y. del N./Bao	-70.79702	19.30579	308	1955-1979	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
67		Sabana Iglesia	Rio Yaque del Norte	Y. del N./Bao	-70.74619	19.31413	230	1967-1980	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
68		El Puente	Rio Yaque del Norte	Y. del N./Dicayagua	-70.76563	19.45218	143	19774-1984	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
69		Guanajuma	Rio Yaque del Norte	del N./Guajuma	-70.74535	19.28108	315	1967-1979	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
70		Jamamú	Rio Yaque del Norte	Y. del N./Jamamu	-70.93344	19.19440	707	1981-1995	Malo	Correo físico: planilla en papel, mensual	si	Instrumentos descalibrados, Irresponsabilidad de los lectores y poca supervisión
71		Jamamú	Rio Yaque del Norte	Y. del N./Jamamu	-70.93344	19.19440	707	1981-1995	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
72		Casa de Máquina	Rio Yaque del Norte	Y. del N./jimenoa	-70.59787	19.09892	602	1981-1987	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
73		Hato Viejo	Rio Yaque del Norte	Y. del N./jimenoa	-70.63414	19.13474	542	1955-1995	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
74		Pinalito	Rio Yaque del Norte	Y. del N/ Jagua	-70.75897	19.29524	283	1955-1979	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
75		Potrero	Rio Yaque del Norte	Y. del N/Amina	-70.96008	19.46969	115	1955-1967	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
76		Paso de la Perra	Rio Yaque del Norte	Yaque del Norte	-70.83966	19.07752	960		Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
77		Pinar Quemado	Rio Yaque del Norte	Yaque del Norte	-70.66513	19.09502	614	1955-2001	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
78		Los Velazquitos	Rio Yaque del Norte	Yaque del Norte	-70.67757	19.21274	375	1981-1994	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
79		Tavera	Rio Yaque del Norte	Yaque del Norte	-70.71535	19.26135	375	1964-1971	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
80		Baitoa	Rio Yaque del Norte	Yaque del Norte	-70.70896	19.32774	215	1957-1966	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
81		Baitoa	Rio Yaque del Norte	Yaque del Norte	-70.70896	19.32774	215	1957-1966	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
82		López	Rio Yaque del Norte	Yaque del Norte	-70.71118	19.33413	190	1961-1966	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
83		López	Rio Yaque del Norte	Yaque del Norte	-70.71118	19.33413	190	1961-1966	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
84		Las Charcas	Rio Yaque del Norte	Yaque del Norte	-70.71396	19.40857	164	1967-1992	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
85		Pte. Viejo Santiago	Rio Yaque del Norte	Yaque del Norte	-70.71174	19.43107	149	1972-1975	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
86		Pte. Viejo Santiago	Rio Yaque del Norte	Yaque del Norte	-70.71174	19.43107	149	1972-1975	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
87		Santiago	Rio Yaque del Norte	Yaque del Norte	-70.70313	19.44135	150	1955-1966	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
88		Santiago	Rio Yaque del Norte	Yaque del Norte	-70.70313	19.44135	150	1955-1966	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
89		Hato Yaque	Rio Yaque del Norte	Yaque del Norte	-70.79619	19.50580	120	1967-1971	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
90		Cañeo	Rio Yaque del Norte	Yaque del Norte	-70.98314	19.54247	80	1967-1989	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
91		Jicome	Rio Yaque del Norte	Yaque del Norte	-71.14731	19.62497	45	1967-1974	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
92		La Caña	Rio Yaque del Norte	Yaque del Norte	-71.52538	19.62247	44	1959-1966	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
93		Dique Juan Gómez	Rio Yaque del Norte	Yaque del Norte	-71.41705	19.68913	23	1975-1994	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
94	42102	El Higuero	Rio Yaque del Norte	Y. del N/ Jagua	-70.83769	19.24795	485	1979-2019	Malo	Telemetría Satelital	si	Estacion Satelital esta cerrada por fuera. Faltan antenas, panel y sensores del rio.
95		Redondo	Rio Yaque del Norte	Y. del N./jimenoa	-70.58389	19.02863	996	1981-1993	Regular	Correo fisico: planilla en papel, mensual	si	Falta la mira #4 y esta tamado de maleza
96	42005	Aguas Caliente	Rio Yaque del Norte	Y. del N./Bao	-70.89814	19.24524		1979-1987	Regular	Correo físico: planilla en papel, mensual	si	Faltan las miras 0 y 1.
97	183101	Pinalito	Rio Yaque del Norte	Y. del N/ Jagua	-70.75897	19.29524		1955-2019	Regular	Correo físico: planilla en papel, mensual	si	Aun se puede considerar de utilidad, porque aun están funcionando y reportando
98	43001	Inoa	Rio Yaque del Norte	Y. del N/Amina	-70.98300	19.35222	348	1967-2017	Regular	Correo físico: planilla en papel, mensual	si	Aun se puede considerar de utilidad, porque aun están funcionando y reportando
99		Paso de la Palma	Rio Yaque del Norte	Y. Del N./Maguaca	-71.51604	19.57455	85	1979-2018	Regular	Correo físico: planilla en papel, mensual	si	Instrumentos descalibrados, Irresponsabilidad de los lectores y poca supervisión
100		Los Pilones	Rio Yaque del Norte	Y. del N. /Baiguaque	-70.77352	19.25579	459	1979-2019	Regular	Correo físico: planilla en papel, mensual	si	Instrumentos descalibrados, Irresponsabilidad de los lectores y poca supervisión
101	42202	El Cerrazo	Rio Yaque del Norte	del N./Guajuma	-70.72710	19.23592	449	1979-2018	Regular	Correo físico: planilla en papel, mensual	si	Faltan las 0 y 3.
102	491005	Palomino	Rio Yaque del Sur	Y. del S./del Medio	-70.96944	18.80512	0	1978-2017	Bueno	Correo físico: planilla en papel, mensual	si	Aun se puede considerar de utilidad, porque aun están funcionando y reportando
103	493006	Paso de Lima	Rio Yaque del Sur	Y. Del S./Río San Juan	-71.29826	19.02472	673	1975-2018	Bueno	Correo físico: planilla en papel, mensual	si	Hay que hacer limpieza del coral. La mira están bien 0-1-2- 3. Punto de mira vieja se esta leyendo
104	493802	Vallejuelo	Rio Yaque del Sur	Y.del S./Los Baos	-71.31721	18.65660	654	1978-2018	Bueno	Correo físico: planilla en papel, mensual	si	La lectora de mira esta enferma
105		La Meseta	Rio Yaque del Sur					1989-2019	Bueno	Correo físico: planilla en papel, mensual		Actualmente están reportando
106		Cucurucho	Rio Yaque del Sur					1986-2019	Bueno	Correo físico: planilla en papel, mensual		Actualmente están reportando
107	493601	Hato Viejo	Rio Yaque del Sur	Y. del S./Arroyo Loro	-71.25538	18.79707	381	1971-2000	Malo	Correo físico: planilla en papel, mensual	si	Existe caseta vieja, La crecida se llevo el linigrafo, No tiene Lector de mira la estación.
108		Rosario	Rio Yaque del Sur	Rio San Juan			381		Malo	Correo físico: planilla en papel, mensual	si	No existe nada, se tomo la coordenada en el punto viejo de la mira. La observadora fue cancelada.
109		Naranjal	Rio Yaque del Sur	Y. del S./Maguana				1964-1967	Malo	Correo físico: planilla en papel, mensual	si	Existe una mira de 3, y el cause del Rio cambio
110		Palo Alto	Rio Yaque del Sur	Río Yaque del Sur	-71.17364	18.29928	0	1967-2011	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
111		Conuquito	Rio Yaque del Sur	Río Yaque del Sur	-71.15148	18.41161	33	1964-1993	Malo	Correo físico: planilla en papel, mensual	si	Instrumentos descalibrados, Irresponsabilidad de los lectores y poca supervisión
112		El Rodeo	Rio Yaque del Sur	Y. Del S./Arroyo Limón	-71.26760	19.00551	0	1982-1994	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
113		La Coja	Rio Yaque del Sur	Y. del S./Guayabal	-70.81868	18.76856	0	1982-2000	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
114		La Guama	Rio Yaque del Sur	Y. del S./Las Cuevas	-70.81296	18.72572	770	1982-2000	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
115		Naranjal	Rio Yaque del Sur	Y. del S./Maguana				1964-1967	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
116	493201	El Cacheo	Rio Yaque del Sur	Y. del S./Mijo				1983-2018	Regular	Correo físico: planilla en papel, mensual		Actualmente están reportando
117		Los Guiros	Rio Yaque del Sur	Río Yaque del Sur	-71.00095	18.51058	0	1985-2018	Regular	Correo físico: planilla en papel, mensual	si	Instrumentos descalibrados, Irresponsabilidad de los lectores y poca supervisión
118		Cañafistol	Rio Yaque del Sur	Rio San Juan			492		Regular	Correo físico: planilla en papel, mensual	si	La mira esta fuera de cauce, se alejo 4 metros. Estación Hidrométrica Satelital.
119		El Cacheo	Rio Yaque del Sur	Y. del S./Arroyo Grande	-71.11519	18.81385	460	1983-2019	Regular	Correo físico: planilla en papel, mensual	si	Hay que cambiar la Mira-Cauce se aleja de la mira.
120		Paso de Lima	Rio Yaque del Sur	Y. Del S./Río San Juan	-71.29826	19.02472	673	1975-2018	Regular	Telemetría Satelital	si	Estación Satelital. Hidrolab no existe, se la llevo la crecida. Existían cuatro mira vieja la cual no sirven. La Nivelación de mira no sirven.
121	100003	Yásica Arriba	Rio Yasica	Yásica/Yásica	-70.59327	19.63365	120	1971-2018	Bueno	Correo fisico: planilla en papel, mensual	si	No tiene observador
122		Palma Herrera	Rio Yasica	Yásica/Partido	-70.41089	19.55469	140	1964-1966	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
123		Sonador	Rio Yasica	Yásica/Sonador	-70.53943	19.59870	133	1956-1984	Malo	Correo físico: planilla en papel, mensual	si	Instrumentos descalibrados, Irresponsabilidad de los lectores y poca supervisión
124	101001	Jamao	Rio Yasica	Yásica/Jamao	-70.44828	19.63171	34	1957-1987	Malo	Telemetría Satelital	si	La ultima crecida del Rio Jamao arrastro las miras.
125	100002	Los Brazos	Rio Yasica	Yásica/Sonador	-70.43089	19.66052	15	1955-1996	Malo	Telemetría Satelital	si	Las miras fueron destruidas con la Tormenta Noel.
126	185502	Santa Ana	Rio Yuna	Yuna/Cenoví	-70.37064	19.25907	82	1982-2019	Bueno	Correo físico: planilla en papel, mensual	si	Actualmente están reportando
127	120002	Yuna	Rio Yuna	Yuna	-70.25278	18.94722	76	2015-2017	Bueno	Telemetría Satelital	Incomp.	Esta todo en buenas condiciones, pero no transmite.
128	181101	Los Arroces	Rio Yuna	Yuna/Masipedro	-70.16638	18.97866	218	1983-2019	Bueno	Telemetría Satelital	si	La observadora de esta estación esta en proceso de pensión, esta muy enferma, por lo que seria conveniente nombrar un observador con urgencia.
129		Blanco (Presa)	Rio Yuna	Rio Blanco	-70.52139	18.88278	465		Malo	Correo físico: planilla en papel, mensual	si	En el muro de la presa de Blanco
130	183001	El Torito	Rio Yuna	Yuna/Blanco	-70.45590	18.75634	760	1982-2004	Malo	Correo físico: planilla en papel, mensual	si	Lleno de maleza dentro y Fuera del corral
131		Los Plátanos	Rio Yuna	Yuna/Maimón	-70.31890	18.78632	240	1958-1965	Malo	Correo físico: planilla en papel, mensual	si	Esta llena de maleza, y no tenemos llave para accesar. Esta transmitiendo.
132		Blanco	Rio Yuna	Yuna/Blanco	-70.52146	18.88300	465	1977-1988	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación. Actualmente no esta en la red de observación
133		Rodeo	Rio Yuna	Yuna/Blanco	-70.56840	18.88439	660	1977-1977	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación. Actualmente no esta en la red de observación
134		Tireo	Rio Yuna	Yuna/Blanco	-70.63285	18.91912	1165	1978-1979	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación. Actualmente no esta en la red de observación
135		Los Guázaros	Rio Yuna	Yuna/Caña Amarga	-70.52646	18.85579	630	1980-1988	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación. Actualmente no esta en la red de observación
136		Tireo	Rio Yuna	Yuna/Tireo	-70.56840	18.88439	660	1980-1995	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del río arrastró las miras
137		El Pino de Yuna	Rio Yuna	Yuna/Yuna	-70.45618	18.79550	582	1982-2006	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida. Tiene un observador nombrado
138		Piedra los Veganos	Rio Yuna	Yuna/Yuna	-70.47173	18.81995	470	1982-2004	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación. Actualmente no esta en la red de observación
139		Hatillo	Rio Yuna	Yuna/Yuna	-70.25283	18.94745		1955-1970	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación. Actualmente no esta en la red de observación
140		Arroyón Arriba	Rio Yuna	Yuna/Arroyón	-70.52827	18.85972	700	1980-1984	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
141		Arroyón Abajo	Rio Yuna	Yuna/Arroyón	-70.52354	18.86066	600	1979-1987	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
142		Jeremías	Rio Yuna	Yuna/Camú	-70.48597	19.22554	70	1957-1961	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
143		Jamo	Rio Yuna	Yuna/Camú	-70.48885	19.22868	65	1973-1978	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
144		Cañabón	Rio Yuna	Yuna/Cañabón	-70.41887	19.05395	118	1964-1966	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
145		Cevicos	Rio Yuna	Yuna/Cevicos	-69.98432	19.00617	78	1963-1965	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
146		Sabana Grande	Rio Yuna	Yuna/Chacuey	-70.06621	19.03554	50	1959-1965	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
147		Jayaco	Rio Yuna	Yuna/Jima	-70.43936	19.00590	161	1955-1965	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
148		Blanco	Rio Yuna	Yuna/Juan de Jesús	-70.36810	18.88717	630	1980-1995	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
149		Naranjal	Rio Yuna	Yuna/Licey	-70.51450	19.32830	96	1964-1987	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
150		Maimón	Rio Yuna	Yuna/Maimón	-70.28490	18.89937	120	1968-2008	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
151		Los Guázaros	Rio Yuna	Yuna/Tireito	-70.52257	18.85886	630	1978-1990	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
152		El Meche	Rio Yuna	Yuna/Toro Flaco	-70.53291	18.84794	700	1980-1986	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
153		Rincón de Yuboa	Rio Yuna	Yuna/Yuboa	-70.36157	18.78331	400	1984-1988	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
154		Los Quemados	Rio Yuna	Yuna/Yuna	-70.47314	18.87916		1962-1984	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
155		Bonao	Rio Yuna	Yuna/Yuna	-70.40762	18.95916	155	1960-1965	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
156		Los Corosos	Rio Yuna	Yuna/Yuna	-70.20089	19.05496	55	1959-1965	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
157		La Jagua	Rio Yuna	Yuna/Yuna	-69.70796	19.14882	13	1997-2004	Malo	Correo físico: planilla en papel, mensual	si	Las crecidas del rio destruyeron la estación, debe ser reconstruida.
158	185003	La Bija	Rio Yuna	Yuna/Camú	-70.12723	19.14948	28	1968-2007	Malo	Telemetría Satelital	si	Esta estación fue vandalizada, solo quedo la caseta. La observadora en proceso de pensión; la señora ya no vive ahí.
159	185001	Bayacanes	Rio Yuna	Yuna/Camú	-70.58153	19.23565	130	1960-1999	Malo	Telemetría Satelital	si	No hay miras. Vandalizada, solo queda la caseta.
160	185501	Cenoví, Los Limones	Rio Yuna	Yuna/Cenoví	-70.34906	19.31966	0	2005-2008	Malo	Telemetría Satelital	si	Las miras fueron arrastradas por una crecida en la Tormenta Noel.
161	184003	El pino de Yuna	Rio Yuna	Yuna/Maimón	-70.29723	18.89472	98	1956-1966	Malo	Telemetría Satelital	si	Esta estación desapareció a raíz de la tormenta Noel.
162	187002	Abadesa	Rio Yuna	Yuna/Payabo	-69.92643	19.01863	51	2015-2017	Malo	Telemetría Satelital	Incomp.	La sección de mira no existen por la crecida del rio, no hay observador, esta pensionado.
163	184002	Los Plátanos	Rio Yuna	Yuna/Yuna	-70.11228	19.12214	39	1971-1982	Malo	Telemetría Satelital	si	Las crecidas del rio destruyeron la estacion. Actualmente no esta en la red de observacion
164	185002	Ranchito	Rio Yuna	Yuna/Camú	-70.38872	19.17508	56	1968-2006	Malo	Telemetría Satelital	si	Tiene un poco de malezas, no funciona.
165	188002	Los Tres Pasos	Rio Yuna	Yuna/Chacuey	-70.07879	18.97640	77	1984-1993	Malo	Telemetría Satelital	si	Las miras no están y no tenemos llave para ver el limnigrafo.
166	188008	La Cabirma	Rio Yuna	Yuna/Maguaca	-70.11522	18.97470	46	1982-2002	Malo	Telemetría Satelital	si	Las miras fueron arrastradas por el rio a raíz de una crecida, hace años. Hay un corral, recién construido, pero no instalaron la telemétrica.
167	183101	Pinalito	Rio Yuna	Yuna/Tireo	-70.28923	18.89662	1165	1981-2019	Malo	Telemetría Satelital	si	No transmite, en mal estado total. Falta mira, además el 50% del corral esta en el suelo.
168		La Boca	Rio Yuna	Yuna/Jaya	-70.30209	19.24578	87	1982-2007	Regular	Correo físico: planilla en papel, mensual	si	Esta estación esta completa, pero no transmite.
169	185603	Jaya	Rio Yuna	Yuna/Jaya	-70.30250	19.24556	61	2015-2017	Regular	Telemetría Satelital	Incomp.	Esta estación esta completa, pero no transmite.
170	185201	Rincón	Rio Yuna	Yuna/Jima	-70.40695	19.10509	80	1957-1976	Regular	Telemetría Satelital	si	Esta estación esta completa, pero no transmite.
171	180003	Villa Riva	Rio Yuna	Yuna/Yuna	-69.90444	19.17440	10	1955-2002	Regular	Telemetría Satelital	si	Quedan 2 miras.
172	180004	El Limón	Rio Yuna	Yuna/Yuna	-69.81921	19.15301	8	1968-2019	Regular	Telemetría Satelital	si	A esta estación le falta 8 metros de miras, por lo que se le dificulta a la observadora hacer bien su trabajo.
173	42301	BAIGUAQUEPILONES LOS	Yaque del Norte	Baiguaque	-70.87945	19.34889	413	1979 - 2018	Bueno	Correo físico: planilla en papel, mensual		Falta la mira 0.
174	42301	BAIGUAQUEPILONES LOS	Yaque del Norte	Baiguaque	-70.87945	19.34889	413	2015-2017	Bueno	Telemetría Satelital	Incomp.	Falta la mira 0.
175		Sabana Alta	Yaque del Sur					1985 - 2019	Bueno	Correo físico: planilla en papel, mensual		Actualmente están reportando

indrhi_hidrometricas_final_sf <- indrhi_hidrometricas_final %>% 
  mutate(excluir_sf = ifelse(is.na(lon_dd_consolidadas) | is.na(lat_dd_consolidadas), 'sí', 'no')) %>% 
  filter(excluir_sf=='no') %>% 
  st_as_sf(coords = c('lon_dd_consolidadas', 'lat_dd_consolidadas'), crs = 4326) %>% 
  select(-excluir_sf) %>%
  mutate(Estado = factor(Estado, levels = c('Bueno','Regular', 'Malo')))

Del total de 175 estaciones hidrométricas, para fines de representación, descartamos 9 que no contaban con coordenadas, o que sus valores eran erróneos. El mapa de estaciones hidrométricas con coordenadas depuradas (166 estaciones), e indicación de sus correspondientes estados, se muestra a continuación.

fpal_estado <- colorFactor(
  palette = c("#E41A1C", "#FFFF33","#4DAF4A"),
  # palette = RColorBrewer::brewer.pal(length(unique(indrhi_hidrometricas_final_sf$Estado)), 'Set1'),
  domain = unique(indrhi_hidrometricas_final_sf$Estado), reverse = T)
leaflet(indrhi_hidrometricas_final_sf) %>%
  addCircleMarkers(
    radius = 5, label = ~`N°`, group = ~ Estado,
    popup = ~ paste0('Codigo N°: ', `N°`, '<br>Nombre: ', `Nombre de la estación`, '<br>Estado: ', Estado),
    color = ~ fpal_estado(Estado),
    stroke = F, fillOpacity = 1
  ) %>%
  addLegend(pal = fpal_estado, values = ~ Estado, opacity = 1,
            title = "Estaciones<br>Hidrométricas<br>INDRHI<br>Estado") %>% 
  addTiles(group = 'OSM') %>%
  addProviderTiles("Esri.NatGeoWorldMap", group="ESRI Mapa") %>%
  addProviderTiles("Esri.WorldImagery", group="ESRI Imagen") %>%
  addProviderTiles("CartoDB.Positron", group= "CartoDB") %>%
  addLayersControl(
    baseGroups = c("ESRI Imagen", "CartoDB", "OSM", "ESRI Mapa"),
    overlayGroups = ~ Estado, position = 'bottomright',
    options = layersControlOptions(collapsed = FALSE)) %>% 
  leaflet_map_view %>% 
  addFullscreenControl()

# Exportar
ruta_gpkg_est_hidrometricas_indrhi <- 'out/con_indicacion_estatus_hidrometricas_indrhi.gpkg'
indrhi_hidrometricas_final_sf %>%
  st_write(ruta_gpkg_est_hidrometricas_indrhi, delete_dsn = TRUE, quiet = T, verbose = F)
ruta_csv_est_hidrometricas_indrhi <- 'out/con_indicacion_estatus_hidrometricas_indrhi.csv'
indrhi_hidrometricas_final %>%
  select(ID = `Código`, Nombre = `Nombre de la estación`, Estado,
         Longitud = lon_dd_consolidadas, Latitud = lat_dd_consolidadas) %>%
  write_csv(ruta_csv_est_hidrometricas_indrhi)

NOTA. Se han generado archivos GeoPackage (out/con_indicacion_estatus_hidrometricas_indrhi.gpkg) y CSV (out/con_indicacion_estatus_hidrometricas_indrhi.csv), de la lista de estaciones hidrométricas de INDRHI, depuradas según el procedimiento explicado en el texto.

Estaciones climáticas

# indrhi <- map(5:6, ~ read_xlsx('fuentes/combinadas/combinadas_v0.9.xlsx', sheet = .x))
indrhi_historico <- read_xlsx('fuentes/indrhi/Listado Red Medicion INDRHI_Historico_24-10-2022_revision_jr.xlsx')
indrhi[[2]] <- indrhi[[2]] %>% rename(Estado = Estatus) %>% 
  filter(!(`N°` == 18 | `N°` == 43)) %>% # Quitar Santana y Jarabacoa, duplicadas
  mutate(Estado = factor(Estado, labels = c('Bueno', 'Malo', 'Regular'))) %>% 
  mutate(Estado = factor(Estado, levels=c('Bueno', 'Regular', 'Malo')))
indrhi_climaticas_resumen_estado <- indrhi[[2]] %>% 
  select(`Código`, `Nombre de la  estación`, Estado) %>% 
  group_by(Estado) %>%
  count()
indrhi_climaticas_resumen_estado %>%
  kable(booktabs=T) %>%
  kable_styling(
    full_width = F,
    latex_options = c("HOLD_position", "scale_down")) %>%
  gsub(' NA ', '', .)

Estado	n
Bueno	17
Regular	12
Malo	32

La fuente empleada para extraer la lista de estaciones climáticas, con sus respectivos estados de operación, fue la misma que la usada para las hidrométricas, es decir, el inventario procedente del PDF Informe Final inventario estaciones hidrometeorológicas INDRHI, Rep. Dom..pdf. La categorización de estaciones según estados indica que, del total de 61 estaciones climáticas, 32 se encontraban en la categoría “Malo”, 12 en estado “Regular” y 17 en la categoría “Bueno”.

La tabla fuente carecía del necesario atributo “coordenadas”, así que, para avanzar en la generación del mapa de estaciones climáticas del INDRHI, nuevamente hubo que manipular los datos (data wrangling). Suplimos la ausencia de coordenadas uniendo la tabla inventario (sin coordenadas) con una lista histórica de estaciones georreferenciadas (mayoritariamente) con coordenadas UTM y/o LAT/LON. Esta última tabla fue aportada por el departamento de hidrología de INDRHI, y corregida puntualmente por nosotros (Listado Red Medicion INDRHI_Historico_24-10-2022_revision_jr.xlsx). Aprovechamos la existencia de campos comunes (códigos, nombre de estación) entre este archivo y la tabla de inventario procedente del PDF. No obstante, tres nuevas dificultades afloraron:

# Unión por código
indrhi[[2]] %>% mutate(Código = as.character(Código)) %>%
  select(CODIGO = Código, ESTACION = `Nombre de la  estación`) %>% 
  inner_join(indrhi_historico %>% select(ORDEN, CODIGO, `X (UTM)`, `Y (UTM)`), by = 'CODIGO') %>% 
  group_by(CODIGO) %>% 
  mutate(dist_x=max(`X (UTM)`) - min(`X (UTM)`), dist_y=max(`Y (UTM)`) - min(`Y (UTM)`)) %>% 
  write_csv('fuentes/indrhi/revision_estaciones_climaticas/para_revision_unidas_por_codigo.csv')
# Resto: no pudieron unirse por código, y se unieron por nombre de estación (búsqueda de texto difusa)
# Probados todos los métodos, "cosine" conservó el máximo posible e hizo las mejores coincidencias
sapply(
  c("osa", "lv", "dl", "hamming", "lcs", "qgram", "cosine", "jaccard", "jw","soundex"),
  function(x)
    indrhi[[2]] %>% mutate(Código = as.character(Código)) %>%
    select(CODIGO = Código, ESTACION = `Nombre de la  estación`) %>% 
    anti_join(indrhi_historico %>% select(ORDEN, CODIGO, `X (UTM)`, `Y (UTM)`), by = 'CODIGO') %>%
    select(-CODIGO) %>% 
    stringdist_left_join(
      y = indrhi_historico %>%
        select(ORDEN, ESTACION, `X (UTM)`, `Y (UTM)`, 'ESTACION',   'NOMBRE_CUE', 'UBICACION'),
      ignore_case = T, max_dist = 0.1, method = x,
      by = 'ESTACION', distance_col='distancia') %>% mutate(metodo=x) %>%
    write_csv(paste0('fuentes/indrhi/revision_estaciones_climaticas/para_revision_unidas_por_nombre_metodo_', x, '.csv')),
  simplify = F)

1. En la tabla de inventario, el campo Código no fue rellenado para todas las estaciones. En este caso, se complementó la unión usando el campo nombre de estación. Usamos el criterio de provincia, cuenca y localidad, para asegurarnos que las estaciones con nombres genéricos (e.g. “CACHEO”) y uso generalizado, fuesen unívocamente identificadas en el par de tablas. Al unir por nombre estación, también se corría el riesgo de que, nombres coincidentes, no fuesen correctamente emparejados al comparar cadenas parecidas pero no exactamente iguales (diferencias de grafías, signos de puntuación, omisiones de artículos). Para ello, evaluamos coincidencia de los campos de nombres de estación en las dos tablas (la de inventario, sin coordenadas, y la histórica, con coordenadas), mediante coincidencia aproximada de cadenas, búsqueda de texto difuso y funciones de distancia de cadenas.

2. En el archivo histórico, muchas estaciones son repeticiones (duplicados, triplicados, …), cada copia con coordenadas iguales y/o ligeramente diferentes entre sí. Normalmente, las distancias entre repeticiones eran de 0 a 100 m, pero también encontramos extremos, como 20 km (ver ejemplo en figura 1), por lo que hubo eligir la instancia más probable según criterio propio, conocimiento del terreno y fuentes cartográficas en entorno SIG (software QGIS).

3. No todas las estaciones disponían de coordenadas en el archivo histórico. Lógicamente, las estaciones sin coordenadas, las excluimos del mapa.

Figure 1: Ejemplos de estaciones repetidas en archivo histórico INDRHI

indrhi_climaticas_depurado <- indrhi[[2]] %>%
               select(ESTACION = `Nombre de la  estación`, Cuenca, `Sub-cuenca`, PROVINCIA,
                      Municipio, `Nivel de Tecnología`, Estado, `Proyecto /Red`, `Observación`) %>%
  left_join(bind_rows(
    read_ods('fuentes/indrhi/revision_estaciones_climaticas/para_revision_unidas_por_codigo.ods') %>% 
      select(ESTACION, ORDEN),
    read_ods('fuentes/indrhi/revision_estaciones_climaticas/para_revision_unidas_por_nombre_metodo_cosine.ods') %>% 
      select(ESTACION = ESTACION.x, ORDEN)),
             by = 'ESTACION') %>% 
  left_join(indrhi_historico %>% select(ORDEN, CODIGO, `X (UTM)`, `Y (UTM)`), by = 'ORDEN')
indrhi_climaticas_depurado_ll <- indrhi_climaticas_depurado %>%
  mutate(excluir_sf = ifelse(is.na(`X (UTM)`) | is.na(`Y (UTM)`), 'sí', 'no')) %>% 
  filter(excluir_sf=='no') %>% select(-excluir_sf) %>% 
  st_as_sf(coords = c('X (UTM)', 'Y (UTM)'), crs = 32619) %>% 
  st_transform(4326) %>% 
  mutate(lon_dd = unlist(map(.$geometry, 1)),
         lat_dd = unlist(map(.$geometry, 2))) %>% 
  st_drop_geometry() %>% 
  select(ESTACION, lon_dd, lat_dd) %>% 
  right_join(indrhi_climaticas_depurado, by = 'ESTACION') %>% 
  relocate(lon_dd:lat_dd, .after = Municipio)
indrhi_climaticas_depurado_ll %>% 
  select(-ORDEN, -`X (UTM)`, -`Y (UTM)`) %>% 
  arrange(ESTACION) %>% 
  kable(booktabs=T) %>%
  kable_styling(
    full_width = F,
    latex_options = c("HOLD_position", "scale_down")) %>%
  gsub(' NA ', '', .)

ESTACION	Cuenca	Sub-cuenca	PROVINCIA	Municipio	lon_dd	lat_dd	Nivel de Tecnología	Estado	Proyecto /Red	Observación	CODIGO
Angelina	Yuna	Yuna	Sánchez Ramírez	Cotui	-70.22229	19.13080	Convencional	Bueno	Red Tradicional Indrhi	Esta estación tiene muchos equipos, aunque funcionan están en mal estado. De la garita no queda nada.	1815
Angostura	Enriquillo	Enriquillo	Independencia	Mella	-71.39485	18.26713	Convencional	Bueno	Red Tradicional Indrhi	La estación esta en malas condiciones	5312
Barraquito	Yuna	Yuna	Duarte	Villa Riva	-69.78891	19.12899	Convencional	Malo	Red Tradicional Indrhi	Esta estación esta en total abandono, no hemos podido contactar al observador.	1814
Bayahibe	Chavón	Chavón	Higuey	Bayahibe	-68.84189	18.40092	Automática	Malo	Red Aquater 1998	Solo existe el corral, no tiene ningún tipo de equipos.	6688
Blanco (Presa)	Yuna	Blanco	Monseñor Noel	Bonao	-70.55812	18.88745	Automática	Regular	Red Banco Mundial 2015 (Sutron)	La estación esta llena de maleza, dentro y fuera del coral, además no transmite.	183201
Boca de Mao	Yaque Del Norte	Yaque Del Norte	Valverde	Mao	-71.05147	19.58830	Convencional	Malo	Red Tradicional Indrhi	Solo funciona el Pluviómetro Mecánico.	406
Catanamatias	Artibonito				-71.40899	19.04968	Convencional	Bueno	Red Tradicional Indrhi		5408
Cenovi en Santa Ana	Yuna	Cenovi en Santa Ana	Duarte	San Francisco de Macorís	-70.34796	19.29467	Convencional/Auto mática	Regular	Red Banco Mundial 2015 (Sutron)	No funciona el tanque evaporímetro	1838
Cítricos Dominicanos	Haina	Guananito	San Cristóbal	Villa Altagracia	-70.24028	18.73897	Automática	Malo	Red Aquater 1998	Solo existe el corral, no tiene ningún tipo de equipos.	6696
Comate	Ozama	Rio Comate	Monte Plata	Bayaguana	-69.61664	18.74568	Automática	Malo	Red Aquater 1998	La estación no existe.	6391
Constanza	Yaque del Sur		La Vega	Constanza	-70.71674	18.91134	Convencional	Bueno	Red Tradicional Indrhi	Es una de las mejores estaciones en la región, de las mas estables.	4902
El Aguacate	Yuna	Yuna	Duarte	Villa Riva	-69.75949	19.17220	Convencional	Malo	Red Tradicional Indrhi	El Sindico de El Aguacate construyo un parque en el área de la estación.	1843
El Arroyaso	Yaque del Norte		La Vega	Jarabacoa	-70.53229	19.04496	Convencional	Malo	Red Tradicional Indrhi	En esta estación solo queda el anemómetro y el pluviómetro. No tiene observador , las observaciones la hacen en el proyecto Ébano Verde.	467
El Limón	Yuna	Yuna	María Trinidad Sánchez	Nagua	-69.81945	19.15301	Automática	Regular	Red Banco Mundial 2015 (Sutron)	El candado de esta estación fue violado, pero todo lo demás esta ahí, tiene un poco de maleza; esta estación no esta transmitiendo.	1803
El Memiso	Rio Ocoa	Rio Ocoa	Azua	Estebania	-70.57159	18.51421	Convencional	Bueno	Red Tradicional Indrhi	Estación libre de vegetación, el corral esta oxidado.	4425
El Naranjito	Yaque Del Norte	Yaque Del Norte	Santiago Rodríguez		-71.34351	19.47442	Convencional	Malo	Red Tradicional Indrhi	La observadora esta pensionada en espera del nombramiento de sustituta.	407
El Seibo	Soco	Rio Soco	El Seibo	El Seibo	-69.04444	18.76550	Convencional	Malo	Red Tradicional Indrhi	No existe ningún equipo en la estación por deterioro y vandalismo, el corral se traslado al patio de la observadora para posible instalación de estación con tecnología Gprs.	3001
Engombe	Haina	Haina	Santo Domingo	Santo Domingo Oeste	-70.00194	18.44941	Convencional	Regular	Red Tradicional Indrhi	Solo funciona el pluviómetro y el anemómetro, el corral en mal estado.	3401
Guarey-Guaigui	Yuna/Camu	Guaigui	La Vega	La Vega			Convencional/Auto mática	Bueno	Red Banco Mundial 2015 (Sutron)	Esta en perfecta condiciones, esta transmitiendo.
Hato Mayor	Higuamo	Higuamo	Hato Mayor	Hierba Buena	-69.25001	18.75023	Automática	Malo	Red Aquater 1998	La estación ya no existe.	3104
Higuerito	Artibonito	Artibonito	Elías Piña	Higuerito	-71.62430	18.99017	Automática	Malo	Red Aquater 1998	Esta estación necesita que se le instale un pluviómetro convencional, ya que el sistema eléctrico no funciona. No tiene el panel solar, ni la tarjeta de memoria; la base del tanque de evaporación esta destruida.	6701
Higuey	Rio Duey	Rio Duey	Higuey	Higuey	-68.69901	18.61992	Automática	Malo	Red Aquater 1998	Ya no existen estructura ni equipos de esta estación.	6392
Isa	Yaque Del Norte	Yaque Del Norte	Santiago	Santiago	-70.74591	19.44608	Convencional	Malo	Red Tradicional Indrhi		404
Jarabacoa	Yaque Del Norte	Yaque Del Norte	La Vega	Jarabacoa	-70.63896	19.13080	Automática	Bueno	Red Banco Mundial 2015 (Sutron)	Esta Transmitiendo.	401
José Contreras	Yuna	Yuna	Hermanas Mirabal	Villa Trina	-70.45006	19.46691	Automática	Regular	Red Banco Mundial 2015 (Sutron)	Esta todo, menos la garita; pero nada esta funcionando.	1816
Juma-Bonao	Yuna	Yuna	Monseñor Noel	Bonao	-70.38618	18.90023	Automática	Bueno	Red Banco Mundial 2015 (Sutron)	La telemétrica esta llena de maleza, el tanque de evaporación no esta funcionando.	1802
Km 11 Carretera Sánchez	Yaque del Sur	Rio San Juan	San Juan de la Maguana	San Juan	-71.15053	18.75813	Convencional	Malo	Red Tradicional Indrhi	Estación en mal estado.	4903
Km 15	Yaque del Sur		Azua	Compostela			Convencional	Malo	Red Tradicional Indrhi	No funciona.
Km 15	Yaque del Sur		Azua	Compostela			Automática	Malo	Red Banco Mundial 2002 (SIAP)
La Antona	Yaque Del Norte	Yaque Del Norte	Montecristy	Guayubin	-71.40288	19.63358	Convencional	Bueno	Red Tradicional Indrhi	La observadora solicita que la estación sea trasladada mas cerca de su residencia.	408
La Guama	Nisibon		Higuey	Nisibon	-68.66119	18.80965	Automática	Malo	Red Aquater 1998	La estación ya no existe.	6691
La Isabela	Bajabonico		Puerto Plata	Isabela	-71.06397	19.83025	Convencional	Bueno	Red Tradicional Indrhi	Observador invalido, lo ayuda su hijo.	602
La Peñita	Este		Hato Mayor		-70.27561	18.45551	Automática	Malo	Red Aquater 1998	Esta estación ya no existe, desconocidos retiraron la estación completa.	380009
Lago Enriquillo	Enriquillo		Independencia	Descubierta	-71.69811	18.55929	Automática	Malo	Red Banco Mundial 2015 (Sutron)	La estación climática tiene todos los equipos en buenas condiciones. No esta funcionando por falta de mantenimiento. El hidrolab instalado para medir el nivel del agua fue sustraído.	5323
Las Barias	Nizao	Nizao	Peravia	Pizarrete			Automática	Regular	Red Banco Mundial 2015 (Sutron)	Estación completa marca Sutron. La base del tanque de evaporación esta deteriorada.
Los Arroyos	Rio Ocoa	Rio Ocoa	San José de Ocoa	Sabana Larga	-70.52916	18.64196	Convencional	Bueno	Red Tradicional Indrhi	Esta libre de vegetación, corral en mal estado.	4426
Los Botados	Yuna	Yuna	Monseñor Noel	Bonao	-70.57674	18.87078	Convencional/Auto mática	Malo	Red Banco Mundial 2015 (Sutron)	Esta llena de maleza, no funcionan, el observador de la estación climática convencional fue cancelado, aun no han nombrado a nadie en su lugar.	1817
Los Jengibres	Rio Boba	Rio Boba	María Trinidad Sánchez	Nagua	-70.04727	19.43913	Convencional	Malo	Red Tradicional Indrhi	La garita colapso y también los termómetros, el tanque de evaporación tiene la base rota, por lo que no esta funcionando.	1501
Los Naranjos	Rio Duey		Higuey	San Rafael del Yuma			Automática	Malo	Red Aquater 1998	La estación ya no existe.
Matagrande	Yaque Del Norte	Yaque Del Norte	Santiago	Sajoma	-70.98758	19.20857	Convencional	Bueno	Red Tradicional Indrhi	Observador incapacitado, la labor la realizar su sobrina.	411
Matayaya	Artibonito	Artibonito	San Juan de la Maguana	Las Matas de Farfán	-71.58844	18.88357	Convencional	Bueno	Red Tradicional Indrhi	La garita no sirve, el termómetro mínima no sirve; la base del tanque de evaporación no sirve.	5401
Medina	Haina	Haina	San Cristóbal	Medina	-70.14450	18.53523	Convencional	Bueno	Red Tradicional Indrhi	Vegetación dentro y fuera, garita en malas condiciones, pluviómetro aceptable, base del pluviómetro deteriorada, al tanque de evaporación le falta el gancho milimétrico, corral en mal estado.	3402
Naranjal	Rio Ocoa	Rio Ocoa	San José de Ocoa	Sabana Larga	-70.47712	18.54298	Convencional	Malo	Red Tradicional Indrhi	No existe estructura ni equipos de la estación, el observador falleció.	4401
Naranjito	Artibonito	Artibonito			-71.49733	19.28413		Bueno	Red Tradicional Indrhi	Falta mantenimiento preventivo	5410
Palo de Caja	Nizao	Nizao	San José de Ocoa		-70.40006	18.53078	Automática	Regular	Red Banco Mundial 2015 (Sutron)	Vegetación dentro del corral, fue limpiada la estación no esta funcionando, el corral tiene ruptura por vandalismo.	3806
Peñón Barahona	Yaque del Sur	Yaque del Sur	Barahona		-71.18722	18.29662	Convencional	Malo	Red Tradicional Indrhi	La estación tiene que ser reparada.	4904
Presa de Montegrande	Yaque del Sur		Barahona				Automática	Malo	Red Banco Mundial 2002 (SIAP)	El pluviómetro esta desnivelado y no tenemos conocimiento de si esta transmitiendo.
Presa de Aguacate	Nizao	Nizao	San Cristóbal	El Cacao			Automática	Regular	Red Banco Mundial 2015 (Sutron)	Estación completa marca Sutron. La estación tiene un poco de maleza en su interior y no se pudo limpiar por no tener las llaves de la estación.
Presa Hatillo	Yuna	Yuna	Sánchez Ramírez	Cotui	-70.25283	18.94745	Automática	Regular	Red Banco Mundial 2015 (Sutron)	Todo en buena condiciones, pero no transmite.	1806
Presa Jiguez	Nizao	Nizao	San José de Ocoa	Jiguey	-70.38255	18.53328	Automática	Regular	Red Banco Mundial 2015 (Sutron)	La estación no tiene el panel solar, le fue robado; la estación no esta funcionando.	380201
Pueblo Nuevo	Ozama		Monte Plata	Majagual	-70.17477	18.43995	Automática	Malo	Red Aquater 1998	La estación ya no existe.	361001
Puerto Escondido	Enriquillo	Las Damas	Independencia	Puerto Escondido	-71.57233	18.32106	Convencional	Malo	Red Tradicional Indrhi	Pluviómetro en mal estado; pero el tanque evaporímetro, el pluviógrafo y el anemómetro están en buen estado.	532001
Quinigua	Yaque Del Norte	Yaque Del Norte	Santiago	Villa González	-70.77368	19.52663	Convencional	Malo	Red Tradicional Indrhi		405
Sabana de la Mar	Este		Hato Mayor	Sabana de la Mar	-69.38908	19.04605	Automática	Regular	Red Aquater 1998	Hay que confirmar su estado.	6693
Sabana Mula	Artibonito	Artibonito			-71.55650	19.03412		Bueno	Red Tradicional Indrhi		5411
Santana	Rio Bani	Rio Bani	Peravia	Bani	-70.21687	18.27898	Automática	Malo	Red Aquater 1998	En esta estación solo existe el corral, no tiene equipos.	6387
Tavera	Yaque Del Norte	Yaque Del Norte	Santiago	Janico	-70.71813	19.28357	Ambas	Malo	Red Tradicional Indrhi		402
Trepada Alta	Rio Yabon	Rio Yanigua	Hato Mayor	El Valle	-69.45930	18.97636	Automática	Malo	Red Aquater 1998	Esta estación no existe.	7320
Triple Ozama	Ozama	Ozama	Monte Plata	Don Juan	-69.94047	18.81030	Automática	Malo	Red Aquater 1998	La estación ya no existe.	6390
Valdesia	Nizao	Nizao	Peravia		-70.28061	18.40856	Automática	Regular	Red Banco Mundial 2015 (Sutron)	Estación completa marca Sutron. El tanque de evaporación tiene la base de madera deteriorada.	3802
Vallejuelo	Yaque del Sur		San Juan de la Maguana	Vallejuelo	-71.33949	18.65774	Convencional	Bueno	Red Tradicional Indrhi	No existe garita y no tiene termómetro.	4915

# sf
indrhi_climaticas_depurado_ll_sf <- indrhi_climaticas_depurado_ll %>% 
  mutate(excluir_sf = ifelse(is.na(lon_dd) | is.na(lat_dd), 'sí', 'no')) %>% 
  filter(excluir_sf=='no') %>% select(-excluir_sf) %>% 
  st_as_sf(coords = c('lon_dd', 'lat_dd'), crs = 4326)

Del total de 61 estaciones hidrométricas, para fines de representación, descartamos 7 que no contaban con coordenadas. El mapa de estaciones climáticas (54 estaciones), e indicación de sus correspondientes estados, se muestra a continuación.

fpal_estado <- colorFactor(
  palette = c("#E41A1C", "#FFFF33","#4DAF4A"),
  # palette = RColorBrewer::brewer.pal(length(unique(indrhi_hidrometricas_final_sf$Estado)), 'Set1'),
  domain = unique(indrhi_climaticas_depurado_ll_sf$Estado), reverse = T)
leaflet(indrhi_climaticas_depurado_ll_sf) %>%
  addCircleMarkers(
    radius = 5, label = ~ESTACION, group = ~ Estado,
    popup = ~ paste0('Codigo N°: ', CODIGO, '<br>Nombre: ', ESTACION, '<br>Estado: ', Estado),
    color = ~ fpal_estado(Estado),
    stroke = F, fillOpacity = 1
  ) %>%
  addLegend(pal = fpal_estado, values = ~ Estado, opacity = 1,
            title = "Estaciones<br>Climáticas<br>INDRHI<br>Estado") %>% 
  addTiles(group = 'OSM') %>%
  addProviderTiles("Esri.NatGeoWorldMap", group="ESRI Mapa") %>%
  addProviderTiles("Esri.WorldImagery", group="ESRI Imagen") %>%
  addProviderTiles("CartoDB.Positron", group= "CartoDB") %>%
  addLayersControl(
    baseGroups = c("ESRI Imagen", "CartoDB", "OSM", "ESRI Mapa"),
    overlayGroups = ~ Estado, position = 'bottomright',
    options = layersControlOptions(collapsed = FALSE)) %>% 
  leaflet_map_view %>% 
  addFullscreenControl()

# Exportar
ruta_gpkg_est_climaticas_indrhi <- 'out/con_indicacion_estatus_climaticas_indrhi.gpkg'
indrhi_climaticas_depurado_ll_sf %>%
  st_write(ruta_gpkg_est_climaticas_indrhi, delete_dsn = T, quiet = T, verbose = F)
ruta_csv_est_climaticas_indrhi <- 'out/con_indicacion_estatus_climaticas_indrhi.csv'
indrhi_climaticas_depurado_ll %>%
  select(Nombre = ESTACION, Estado,
         Longitud = lon_dd, Latitud = lat_dd) %>%
  write_csv(ruta_csv_est_climaticas_indrhi)

NOTA. Se han generado archivos GeoPackage (out/con_indicacion_estatus_climaticas_indrhi.gpkg) y CSV (out/con_indicacion_estatus_climaticas_indrhi.csv), de la lista de estaciones climáticas de INDRHI, depuradas según el procedimiento explicado en el texto.

Guakía Ambiente

Sin cambios respecto de la primera evaluación realizada en noviembre de 2022 (ver cuaderno consolidacion-lista-estaciones.md.

guakia_sf <- invisible(st_read('out/consolidado_estaciones_sf.gpkg', optional = T, quiet = T)) %>% 
  st_set_geometry('geom') %>% 
  filter(grepl('guakia', fuente, ignore.case = T)) %>% 
  select(Nombre = idOK) %>%
  mutate(Estado = 'activa o bueno', entidad = 'Guakía Ambiente',
         tipo = 'meteoclimática', propiedad = 'privada',
         Longtiud = st_coordinates(geom)[,1], Latitud = st_coordinates(geom)[,2])
guakia_sf %>% st_write('out/con_indicacion_estatus_guakia_climaticas.gpkg',
                       delete_dsn = T, quiet = T, verbose = F)

REDDOM

Sin cambios respecto de la primera evaluación realizada en noviembre de 2022 (ver cuaderno consolidacion-lista-estaciones.md.

reddom_sf <- invisible(st_read('out/consolidado_estaciones_sf.gpkg', optional = T, quiet = T)) %>% 
  st_set_geometry('geom') %>% 
  filter(grepl('reddom', fuente, ignore.case = T)) %>% 
  select(Nombre = idOK) %>%
  mutate(Estado = 'activa o bueno', entidad = 'Fundación REDDOM',
         tipo = 'meteoclimática', propiedad = 'privada',
         Longtiud = st_coordinates(geom)[,1], Latitud = st_coordinates(geom)[,2])
reddom_sf %>% st_write('out/con_indicacion_estatus_reddom_climaticas.gpkg',
                       delete_dsn = T, quiet = T, verbose = F)

INTEC

Sin cambios respecto de la primera evaluación realizada en noviembre de 2022 (ver cuaderno consolidacion-lista-estaciones.md.

intec_sf <- invisible(st_read('out/consolidado_estaciones_sf.gpkg', optional = T, quiet = T)) %>% 
  st_set_geometry('geom') %>% 
  filter(grepl('intec', fuente, ignore.case = T)) %>% 
  select(Nombre = idOK) %>%
  mutate(Estado = 'activa o bueno', entidad = 'INTEC',
         tipo = 'meteoclimática', propiedad = 'privada',
         Longtiud = st_coordinates(geom)[,1], Latitud = st_coordinates(geom)[,2])
intec_sf %>% st_write('out/con_indicacion_estatus_intec_climaticas.gpkg',
                      delete_dsn = T, quiet = T, verbose = F)

PUCMM

Obtuvimos coordenadas por medio de correo electrónico, fechado 19 de diciembre de 2022.

pucmm_df <- data.frame(lon = -70.682416, lat = 19.4409224)
pucmm_sf <- st_as_sf(pucmm_df, coords = c("lon","lat"), crs = 4326) %>% 
  st_set_geometry('geom') %>% 
  mutate(Nombre = 'Campus PUCCM', Estado = 'activa o bueno', entidad = 'INTEC',
         tipo = 'meteoclimática', propiedad = 'privada',
         Longtiud = st_coordinates(geom)[,1], Latitud = st_coordinates(geom)[,2])
pucmm_sf %>% st_write('out/con_indicacion_estatus_pucmm_climaticas.gpkg',
                      delete_dsn = T, quiet = T, verbose = F)

CAEI

Obtuvimos coordenadas de una estación meteoclimática durante reunión virtual en diciembre de 2022 y, posteriormente, coordenadas de una red de pluviómetros por correo electrónico.

caei_pluv_sf <- invisible(st_read('fuentes/caei/PluviometrosCAEI/PluviometrosCAEI.gpkg', optional = T, quiet = T)) %>%
  st_geometry() %>% st_geometry() %>% st_as_sf() %>% st_set_geometry('geom') %>% 
  mutate(tipo = 'pluviométrica') %>% st_zm()
caei_clim_df <- data.frame(lon = -(69+31/60+29.59/3600), lat = 18+35/60+37.83/3600)
caei_clim_sf <- st_as_sf(caei_clim_df, coords = c("lon","lat"), crs = 4326) %>%
  st_geometry() %>% st_as_sf() %>% st_set_geometry('geom') %>% 
  mutate(tipo = 'meteoclimática')
caei_sf <-  bind_rows(caei_pluv_sf, caei_clim_sf) %>% 
  mutate(Estado = 'activa o bueno', entidad = 'CAEI', propiedad = 'privada',
         Longtiud = st_coordinates(geom)[,1], Latitud = st_coordinates(geom)[,2])
caei_sf %>% st_write('out/con_indicacion_estatus_caei_climaticas_pluviometricas.gpkg',
                     delete_dsn = T, quiet = T, verbose = F)

Universidad ISA

isa_df <- data.frame(lon = -70.7482792, lat = 19.4431974)
isa_sf <- st_as_sf(isa_df, coords = c("lon","lat"), crs = 4326) %>% 
  st_set_geometry('geom') %>% 
  mutate(Nombre = 'Estación Japonesa', Estado = 'activa o bueno', entidad = 'Universidad ISA',
         tipo = 'meteoclimática', propiedad = 'privada',
         Longtiud = st_coordinates(geom)[,1], Latitud = st_coordinates(geom)[,2])
isa_sf %>% st_write('out/con_indicacion_estatus_isa_climaticas.gpkg',
                    delete_dsn = T, quiet = T, verbose = F)

Servicio Geológico Nacional (SGN)

sgn_df <- data.frame(x = c(417006, 392375), y = c(2053599, 2042564))
sgn_sf <- st_as_sf(sgn_df, coords = c("x","y"), crs = 32619) %>% 
  st_transform(4326) %>% st_set_geometry('geom') %>% 
  mutate(Nombre = c('Río Ozama en barrera de salinidad', 'Río Haina en Manoguayabo'),
         Estado = 'activa o bueno', entidad = 'SGN',
         tipo = 'hidrométrica', propiedad = 'pública',
         Longtiud = st_coordinates(geom)[,1], Latitud = st_coordinates(geom)[,2])
sgn_sf %>% st_write('out/con_indicacion_estatus_sgn_hidrometricas.gpkg',
                    delete_dsn = T, quiet = T, verbose = F)

EGEHAINA

Aunque se recibió documentación de las estaciones, no se recibieron coordenadas.

JAD

Aunque se recibió documentación de las estaciones, no se recibieron coordenadas.

Consolidado

indrhi_p_consolidado_sf <- bind_rows(
  indrhi_climaticas_depurado_ll_sf %>%
    st_set_geometry('geom') %>% 
    select(Nombre = ESTACION, Estado) %>%
    mutate(tipo = 'meteoclimática'),
  indrhi_hidrometricas_final_sf %>%
    st_set_geometry('geom') %>% 
    select(Nombre = `Nombre de la estación`, Estado) %>% 
    mutate(tipo = 'hidrométrica')) %>% 
  mutate(entidad = 'INDRHI', propiedad = 'pública', Estado = as.character(Estado),
         Longtiud = st_coordinates(geom)[,1], Latitud = st_coordinates(geom)[,2]) %>% 
  mutate(Estado = case_when(
      Estado == 'Bueno' ~ 'activa o bueno',
      Estado == 'Regular' ~ 'regular',
      Estado == 'Malo' ~ 'inactiva o no reportada',
      TRUE ~ Estado))
onamet_p_consolidado_sf <- lista_con_periodo_2022_coord_todas_sf %>% 
  st_set_geometry('geom') %>% st_as_sf() %>% 
  select(Nombre = nombre_onamet, Estado = estado) %>% 
  mutate(tipo = 'meteoclimática', entidad = 'ONAMET', propiedad = 'pública', Estado = as.character(Estado),
         Longtiud = st_coordinates(geom)[,1], Latitud = st_coordinates(geom)[,2]) %>% 
  mutate(Estado = case_when(
      Estado == 'activa' ~ 'activa o bueno',
      TRUE ~ Estado))
consolidado_sf <- bind_rows(
  guakia_sf, intec_sf, pucmm_sf, reddom_sf, caei_sf, isa_sf, sgn_sf,
  indrhi_p_consolidado_sf, onamet_p_consolidado_sf) %>% 
  mutate(Estado = factor(Estado, levels = c('activa o bueno','regular', 'inactiva o no reportada')))
consolidado_sf %>%
  st_write('out/con_indicacion_estatus_consolidado.gpkg',
           delete_dsn = T, quiet = T, verbose = F)
consolidado_sf %>% st_drop_geometry %>% write_csv('out/con_indicacion_estatus_consolidado.csv')

Mapa de estaciones climáticas según estado

fpal_estado <- colorFactor(
  palette = c("#E41A1C", "#FFFF33","#4DAF4A"),
  domain = unique(consolidado_sf$Estado), reverse = T)
leaflet(consolidado_sf %>% filter(tipo == 'meteoclimática')) %>%
  addCircleMarkers(
    radius = 5, group = ~ Estado,
    popup = ~ paste0('Entidad: ', entidad, '<br>Nombre: ', Nombre, '<br>Estado: ', Estado),
    color = ~ fpal_estado(Estado),
    stroke = F, fillOpacity = 1
  ) %>%
  addLegend(pal = fpal_estado, values = ~ Estado, opacity = 1,
            title = "Estaciones<br>Climáticas<br>Estado") %>% 
  addTiles(group = 'OSM') %>%
  addProviderTiles("Esri.NatGeoWorldMap", group="ESRI Mapa") %>%
  addProviderTiles("Esri.WorldImagery", group="ESRI Imagen") %>%
  addProviderTiles("CartoDB.Positron", group= "CartoDB") %>%
  addLayersControl(
    baseGroups = c("ESRI Imagen", "CartoDB", "OSM", "ESRI Mapa"),
    overlayGroups = ~ Estado, position = 'bottomright',
    options = layersControlOptions(collapsed = FALSE)) %>% 
  leaflet_map_view %>% 
  addFullscreenControl()

Mapa de estaciones hidrométricas según estado

leaflet(consolidado_sf %>% filter(tipo == 'hidrométrica')) %>%
  addCircleMarkers(
    radius = 5, group = ~ Estado,
    popup = ~ paste0('Entidad: ', entidad, '<br>Nombre: ', Nombre, '<br>Estado: ', Estado),
    color = ~ fpal_estado(Estado),
    stroke = F, fillOpacity = 1
  ) %>%
  addLegend(pal = fpal_estado, values = ~ Estado, opacity = 1,
            title = "Estaciones<br>Hidrométricas<br>Estado") %>% 
  addTiles(group = 'OSM') %>%
  addProviderTiles("Esri.NatGeoWorldMap", group="ESRI Mapa") %>%
  addProviderTiles("Esri.WorldImagery", group="ESRI Imagen") %>%
  addProviderTiles("CartoDB.Positron", group= "CartoDB") %>%
  addLayersControl(
    baseGroups = c("ESRI Imagen", "CartoDB", "OSM", "ESRI Mapa"),
    overlayGroups = ~ Estado, position = 'bottomright',
    options = layersControlOptions(collapsed = FALSE)) %>% 
  leaflet_map_view %>% 
  addFullscreenControl()

Mapa de estaciones climáticas según propiedad

fpal_propiedad <- colorFactor(
  palette = c("#f1a340", "#998ec3"),
  domain = unique(consolidado_sf$propiedad))
leaflet(consolidado_sf %>% filter(tipo == 'meteoclimática')) %>%
  addCircleMarkers(
    radius = 5, group = ~ propiedad,
    popup = ~ paste0('Entidad: ', entidad, '<br>Nombre: ', Nombre, '<br>Estado: ', Estado),
    color = ~ fpal_propiedad(propiedad),
    stroke = F, fillOpacity = 1
  ) %>%
  addLegend(pal = fpal_propiedad, values = ~ propiedad, opacity = 1,
            title = "Estaciones<br>Climáticas<br>Propiedad") %>% 
  addTiles(group = 'OSM') %>%
  addProviderTiles("Esri.NatGeoWorldMap", group="ESRI Mapa") %>%
  addProviderTiles("Esri.WorldImagery", group="ESRI Imagen") %>%
  addProviderTiles("CartoDB.Positron", group= "CartoDB") %>%
  addLayersControl(
    baseGroups = c("ESRI Imagen", "CartoDB", "OSM", "ESRI Mapa"),
    overlayGroups = ~ propiedad, position = 'bottomright',
    options = layersControlOptions(collapsed = FALSE)) %>% 
  leaflet_map_view %>% 
  addFullscreenControl()

Mapa de estaciones hidrométricas según propiedad (Yaque vendría aquí)

leaflet(consolidado_sf %>% filter(tipo == 'hidrométrica')) %>%
  addCircleMarkers(
    radius = 5, group = ~ propiedad,
    popup = ~ paste0('Entidad: ', entidad, '<br>Nombre: ', Nombre, '<br>Estado: ', Estado),
    color = ~ fpal_propiedad(propiedad),
    stroke = F, fillOpacity = 1
  ) %>%
  addLegend(pal = fpal_propiedad, values = ~ propiedad, opacity = 1,
            title = "Estaciones<br>Climáticas<br>Propiedad") %>% 
  addTiles(group = 'OSM') %>%
  addProviderTiles("Esri.NatGeoWorldMap", group="ESRI Mapa") %>%
  addProviderTiles("Esri.WorldImagery", group="ESRI Imagen") %>%
  addProviderTiles("CartoDB.Positron", group= "CartoDB") %>%
  addLayersControl(
    baseGroups = c("ESRI Imagen", "CartoDB", "OSM", "ESRI Mapa"),
    overlayGroups = ~ propiedad, position = 'bottomright',
    options = layersControlOptions(collapsed = FALSE)) %>% 
  leaflet_map_view %>% 
  addFullscreenControl()

Mapa de estaciones climáticas según entidad

fpal_entidad <- colorFactor(
  palette = RColorBrewer::brewer.pal(length(unique(consolidado_sf$entidad)), 'Set1'),
  domain = unique(consolidado_sf$entidad))
leaflet(consolidado_sf %>% filter(tipo == 'meteoclimática')) %>%
  addCircleMarkers(
    radius = 5, group = ~ entidad,
    popup = ~ paste0('Entidad: ', entidad, '<br>Nombre: ', Nombre, '<br>Estado: ', Estado),
    color = ~ fpal_entidad(entidad),
    stroke = F, fillOpacity = 1
  ) %>%
  addLegend(pal = fpal_entidad, values = ~ entidad, opacity = 1,
            title = "Estaciones<br>Climáticas<br>Entidad") %>% 
  addTiles(group = 'OSM') %>%
  addProviderTiles("Esri.NatGeoWorldMap", group="ESRI Mapa") %>%
  addProviderTiles("Esri.WorldImagery", group="ESRI Imagen") %>%
  addProviderTiles("CartoDB.Positron", group= "CartoDB") %>%
  addLayersControl(
    baseGroups = c("ESRI Imagen", "CartoDB", "OSM", "ESRI Mapa"),
    overlayGroups = ~ entidad, position = 'bottomright',
    options = layersControlOptions(collapsed = FALSE)) %>% 
  leaflet_map_view %>% 
  addFullscreenControl()

Mapa de estaciones hidrométricas según entidad

leaflet(consolidado_sf %>% filter(tipo == 'hidrométrica')) %>%
  addCircleMarkers(
    radius = 5, group = ~ entidad,
    popup = ~ paste0('Entidad: ', entidad, '<br>Nombre: ', Nombre, '<br>Estado: ', Estado),
    color = ~ fpal_entidad(entidad),
    stroke = F, fillOpacity = 1
  ) %>%
  addLegend(pal = fpal_entidad, values = ~ entidad, opacity = 1,
            title = "Estaciones<br>Hidrométricas<br>Entidad") %>% 
  addTiles(group = 'OSM') %>%
  addProviderTiles("Esri.NatGeoWorldMap", group="ESRI Mapa") %>%
  addProviderTiles("Esri.WorldImagery", group="ESRI Imagen") %>%
  addProviderTiles("CartoDB.Positron", group= "CartoDB") %>%
  addLayersControl(
    baseGroups = c("ESRI Imagen", "CartoDB", "OSM", "ESRI Mapa"),
    overlayGroups = ~ entidad, position = 'bottomright',
    options = layersControlOptions(collapsed = FALSE)) %>% 
  leaflet_map_view %>% 
  addFullscreenControl()

Mapas individuales de estaciones climáticas por entidades según estado

ruta_ez_gh <- 'https://raw.githubusercontent.com/geofis/zonal-statistics/'
ez_ver <- 'd7f79365168e688f0d78f521e53fbf2da19244ef/'
pais_url <- paste0(ruta_ez_gh, ez_ver, 'inst/extdata/dr.gpkg')
pais <- invisible(st_read(pais_url, optional = T, layer = 'pais', quiet = T))
st_geometry(pais) <- "geometry"
map(unique(consolidado_sf %>% filter(tipo == 'meteoclimática') %>% pull(entidad)),
    function(x) {
      consolidado_sf %>% filter(entidad == x & tipo == 'meteoclimática') %>% 
      ggplot +
        geom_sf(data = pais, fill = 'transparent', color = 'grey50') +
        geom_sf(alpha = 0.8, aes(fill = Estado), shape = 21, size = 1.5) +
        scale_fill_manual(values = c('activa o bueno'='green', 'regular' = 'yellow', 'inactiva o no reportada' = 'red')) +
        labs(title = paste0('Estaciones climáticas de ', x)) +
       theme_bw() + 
       ggspatial::annotation_scale(style = 'ticks') +
       theme(legend.title = element_blank()) 
    })

## [[1]]

## 
## [[2]]

## 
## [[3]]

## 
## [[4]]

## 
## [[5]]

## 
## [[6]]

## 
## [[7]]

Mapa individuales de estaciones pluviométricas por entidades según estado

map(unique(consolidado_sf %>% filter(tipo == 'pluviométrica') %>% pull(entidad)),
    function(x) {
      consolidado_sf %>% filter(entidad == x & tipo == 'pluviométrica') %>% 
      ggplot +
        geom_sf(data = pais, fill = 'transparent', color = 'grey50') +
        geom_sf(alpha = 0.8, aes(fill = Estado), shape = 21, size = 1.5) +
        scale_fill_manual(values = c('activa o bueno'='green', 'regular' = 'yellow', 'inactiva o no reportada' = 'red')) +
        labs(title = paste0('Estaciones pluviométricas de ', x)) +
       theme_bw() + 
       ggspatial::annotation_scale(style = 'ticks') +
       theme(legend.title = element_blank()) 
    })

## [[1]]

Mapas individuales de estaciones hidrométricas por entidades según estado

map(unique(consolidado_sf %>% filter(tipo == 'hidrométrica') %>% pull(entidad)),
    function(x) {
      consolidado_sf %>% filter(entidad == x & tipo == 'hidrométrica') %>% 
      ggplot +
        geom_sf(data = pais, fill = 'transparent', color = 'grey50') +
        geom_sf(alpha = 0.8, aes(fill = Estado), shape = 21, size = 1.5) +
        scale_fill_manual(values = c('activa o bueno'='green', 'regular' = 'yellow', 'inactiva o no reportada' = 'red')) +
        labs(title = paste0('Estaciones hidrométricas de ', x)) +
       theme_bw() + 
       ggspatial::annotation_scale(style = 'ticks') +
       theme(legend.title = element_blank()) 
    })

## [[1]]

## 
## [[2]]

Impresión de tabla resumen

consolidado_sf %>% st_drop_geometry %>%
  select(entidad, tipo, Estado) %>% 
  group_by_all() %>% 
  count() %>% 
  pivot_wider(names_from = tipo, values_from = n, values_fill = 0) %>% 
  adorn_totals(c('row', 'col')) %>% 
  kable(booktabs=T) %>%
  kable_styling(latex_options = c("HOLD_position", "scale_down")) %>%
  gsub(' NA ', '', .)

entidad	Estado	meteoclimática	pluviométrica	hidrométrica	Total
CAEI	activa o bueno	1	95	0	96
Fundación REDDOM	activa o bueno	32	0	0	32
Guakía Ambiente	activa o bueno	1	0	0	1
INDRHI	activa o bueno	16	0	26	42
INDRHI	regular	10	0	22	32
INDRHI	inactiva o no reportada	28	0	118	146
INTEC	activa o bueno	6	0	0	6
ONAMET	activa o bueno	36	0	0	36
ONAMET	inactiva o no reportada	51	0	0	51
SGN	activa o bueno	0	0	2	2
Universidad ISA	activa o bueno	1	0	0	1
Total		182	95	168	445

Impresión de tabla exhaustiva

consolidado_sf %>% st_drop_geometry %>%
  arrange(entidad) %>% 
  kable(booktabs=T) %>%
  kable_styling(latex_options = c("HOLD_position", "scale_down")) %>%
  gsub(' NA ', '', .)

Nombre	Estado	entidad	tipo	propiedad	Longtiud	Latitud
	activa o bueno	CAEI	pluviométrica	privada	-69.21889	18.55294
	activa o bueno	CAEI	pluviométrica	privada	-69.18947	18.55230
	activa o bueno	CAEI	pluviométrica	privada	-69.19707	18.57829
	activa o bueno	CAEI	pluviométrica	privada	-69.17892	18.59604
	activa o bueno	CAEI	pluviométrica	privada	-69.17080	18.58076
	activa o bueno	CAEI	pluviométrica	privada	-69.32054	18.52682
	activa o bueno	CAEI	pluviométrica	privada	-69.32006	18.50352
	activa o bueno	CAEI	pluviométrica	privada	-69.53086	18.68998
	activa o bueno	CAEI	pluviométrica	privada	-69.59511	18.69934
	activa o bueno	CAEI	pluviométrica	privada	-69.45967	18.49031
	activa o bueno	CAEI	pluviométrica	privada	-69.58603	18.69027
	activa o bueno	CAEI	pluviométrica	privada	-69.57901	18.67273
	activa o bueno	CAEI	pluviométrica	privada	-69.61238	18.62924
	activa o bueno	CAEI	pluviométrica	privada	-69.54196	18.65939
	activa o bueno	CAEI	pluviométrica	privada	-69.55103	18.63020
	activa o bueno	CAEI	pluviométrica	privada	-69.50963	18.58394
	activa o bueno	CAEI	pluviométrica	privada	-69.48853	18.49636
	activa o bueno	CAEI	pluviométrica	privada	-69.33381	18.46710
	activa o bueno	CAEI	pluviométrica	privada	-69.33716	18.46023
	activa o bueno	CAEI	pluviométrica	privada	-69.36816	18.48667
	activa o bueno	CAEI	pluviométrica	privada	-69.23283	18.51922
	activa o bueno	CAEI	pluviométrica	privada	-69.55671	18.60408
	activa o bueno	CAEI	pluviométrica	privada	-69.48490	18.57442
	activa o bueno	CAEI	pluviométrica	privada	-69.49747	18.50183
	activa o bueno	CAEI	pluviométrica	privada	-69.23041	18.53580
	activa o bueno	CAEI	pluviométrica	privada	-69.22398	18.65683
	activa o bueno	CAEI	pluviométrica	privada	-69.62645	18.68470
	activa o bueno	CAEI	pluviométrica	privada	-69.59897	18.67025
	activa o bueno	CAEI	pluviométrica	privada	-69.18396	18.60210
	activa o bueno	CAEI	pluviométrica	privada	-69.32054	18.52682
	activa o bueno	CAEI	pluviométrica	privada	-69.56018	18.69046
	activa o bueno	CAEI	pluviométrica	privada	-69.27009	18.59312
	activa o bueno	CAEI	pluviométrica	privada	-69.46264	18.51042
	activa o bueno	CAEI	pluviométrica	privada	-69.43877	18.51227
	activa o bueno	CAEI	pluviométrica	privada	-69.39803	18.51268
	activa o bueno	CAEI	pluviométrica	privada	-69.53803	18.59163
	activa o bueno	CAEI	pluviométrica	privada	-69.52501	18.59315
	activa o bueno	CAEI	pluviométrica	privada	-69.22686	18.53883
	activa o bueno	CAEI	pluviométrica	privada	-69.53132	18.66737
	activa o bueno	CAEI	pluviométrica	privada	-69.51958	18.63589
	activa o bueno	CAEI	pluviométrica	privada	-69.52051	18.59256
	activa o bueno	CAEI	pluviométrica	privada	-69.27129	18.52132
	activa o bueno	CAEI	pluviométrica	privada	-69.55889	18.62121
	activa o bueno	CAEI	pluviométrica	privada	-69.56821	18.66182
	activa o bueno	CAEI	pluviométrica	privada	-69.60786	18.58096
	activa o bueno	CAEI	pluviométrica	privada	-69.23799	18.54002
	activa o bueno	CAEI	pluviométrica	privada	-70.17930	18.32100
	activa o bueno	CAEI	pluviométrica	privada	-70.20280	18.30830
	activa o bueno	CAEI	pluviométrica	privada	-70.19860	18.33960
	activa o bueno	CAEI	pluviométrica	privada	-70.22170	18.32930
	activa o bueno	CAEI	pluviométrica	privada	-70.22770	18.27830
	activa o bueno	CAEI	pluviométrica	privada	-70.24450	18.25720
	activa o bueno	CAEI	pluviométrica	privada	-70.24698	18.24269
	activa o bueno	CAEI	pluviométrica	privada	-70.18782	18.35632
	activa o bueno	CAEI	pluviométrica	privada	-70.16364	18.33952
	activa o bueno	CAEI	pluviométrica	privada	-69.62263	18.64301
	activa o bueno	CAEI	pluviométrica	privada	-69.41818	18.54273
	activa o bueno	CAEI	pluviométrica	privada	-69.48588	18.55482
	activa o bueno	CAEI	pluviométrica	privada	-69.42933	18.55270
	activa o bueno	CAEI	pluviométrica	privada	-69.42797	18.57315
	activa o bueno	CAEI	pluviométrica	privada	-69.40927	18.57943
	activa o bueno	CAEI	pluviométrica	privada	-69.43684	18.64077
	activa o bueno	CAEI	pluviométrica	privada	-69.43904	18.67101
	activa o bueno	CAEI	pluviométrica	privada	-69.42010	18.61537
	activa o bueno	CAEI	pluviométrica	privada	-69.44632	18.59368
	activa o bueno	CAEI	pluviométrica	privada	-69.45662	18.69537
	activa o bueno	CAEI	pluviométrica	privada	-69.44525	18.73182
	activa o bueno	CAEI	pluviométrica	privada	-69.43675	18.70216
	activa o bueno	CAEI	pluviométrica	privada	-69.41941	18.68300
	activa o bueno	CAEI	pluviométrica	privada	-69.39473	18.69982
	activa o bueno	CAEI	pluviométrica	privada	-69.40992	18.71576
	activa o bueno	CAEI	pluviométrica	privada	-69.42773	18.74102
	activa o bueno	CAEI	pluviométrica	privada	-69.39123	18.66266
	activa o bueno	CAEI	pluviométrica	privada	-69.36895	18.62069
	activa o bueno	CAEI	pluviométrica	privada	-69.37292	18.68376
	activa o bueno	CAEI	pluviométrica	privada	-69.33940	18.57404
	activa o bueno	CAEI	pluviométrica	privada	-69.31271	18.58698
	activa o bueno	CAEI	pluviométrica	privada	-69.31144	18.54462
	activa o bueno	CAEI	pluviométrica	privada	-69.31994	18.61152
	activa o bueno	CAEI	pluviométrica	privada	-69.31324	18.63391
	activa o bueno	CAEI	pluviométrica	privada	-69.32827	18.59361
	activa o bueno	CAEI	pluviométrica	privada	-69.32952	18.64137
	activa o bueno	CAEI	pluviométrica	privada	-69.21798	18.57663
	activa o bueno	CAEI	pluviométrica	privada	-69.18480	18.67372
	activa o bueno	CAEI	pluviométrica	privada	-69.20693	18.51756
	activa o bueno	CAEI	pluviométrica	privada	-69.24272	18.51155
	activa o bueno	CAEI	pluviométrica	privada	-69.40861	18.55778
	activa o bueno	CAEI	pluviométrica	privada	-69.50196	18.50637
	activa o bueno	CAEI	pluviométrica	privada	-69.51059	18.48191
	activa o bueno	CAEI	pluviométrica	privada	-69.46276	18.45942
	activa o bueno	CAEI	pluviométrica	privada	-69.42882	18.59325
	activa o bueno	CAEI	pluviométrica	privada	-69.38022	18.46860
	activa o bueno	CAEI	pluviométrica	privada	-69.40572	18.52338
	activa o bueno	CAEI	pluviométrica	privada	-69.38824	18.51864
	activa o bueno	CAEI	pluviométrica	privada	-69.52489	18.59384
	activa o bueno	CAEI	meteoclimática	privada	-69.52489	18.59384
Los_Tocones_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.34635	19.50247
Chacuey_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.56480	19.50264
Agua_de_Luis_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.21466	19.75922
Juliana_Jaramillo_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.61470	19.78727
Banamiel_Oficina_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.19435	19.66831
Esnamarena_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.60576	19.12550
Unisa_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.74848	19.44292
Agrofrontera_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.38353	19.60421
Los_Montones_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.92510	19.28876
La_Guama_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.50160	18.71723
Plan_Yaque_Manabao_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.79701	19.07087
Ekoban1_Amina_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.98673	19.53784
Banelino_Mao	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.09064	19.55873
Montecristi_Oficina_Banelino	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.64626	19.84958
Amina_Banelino	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.96859	19.54510
Judea_Nueva_Banamiel	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.65865	19.70878
Hatillo_Palma_Banelino	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.21455	19.66124
Hato_al_Medio_Banelino	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.31737	19.67714
Bananera_los_Rios	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.19988	19.63595
Novillero_Clay_REDDOM	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.11765	19.86363
Guananico_REDDOM	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.93845	19.73703
Puerto Plata_REDDOM	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.64679	19.70652
Jumunuco_REDDOM	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.74990	19.09705
La_Pita_REDDOM	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.63360	19.00903
Ucateci_LaVega_REDDOM	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.51910	19.23312
Las_Matas_de_Farfan_REDDOM	activa o bueno	Fundación REDDOM	meteoclimática	privada	-71.54292	18.87926
La_Cumbre_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.62268	19.54729
Santiago_Apedi_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.70682	19.47227
Cevicos_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-69.97047	19.00639
Cimpa_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.84041	19.52976
Conacado_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.25831	19.04199
Sur_Futuro_Reddom	activa o bueno	Fundación REDDOM	meteoclimática	privada	-70.94145	18.73675
Vuelta Larga	activa o bueno	Guakía Ambiente	meteoclimática	privada	-69.99350	19.29050
Catanamatias	activa o bueno	INDRHI	meteoclimática	pública	-71.40899	19.04968
Matayaya	activa o bueno	INDRHI	meteoclimática	pública	-71.58844	18.88357
Naranjito	activa o bueno	INDRHI	meteoclimática	pública	-71.49733	19.28413
Sabana Mula	activa o bueno	INDRHI	meteoclimática	pública	-71.55650	19.03412
Higuerito	inactiva o no reportada	INDRHI	meteoclimática	pública	-71.62430	18.99017
La Isabela	activa o bueno	INDRHI	meteoclimática	pública	-71.06397	19.83025
Bayahibe	inactiva o no reportada	INDRHI	meteoclimática	pública	-68.84189	18.40092
Angostura	activa o bueno	INDRHI	meteoclimática	pública	-71.39485	18.26713
Lago Enriquillo	inactiva o no reportada	INDRHI	meteoclimática	pública	-71.69811	18.55929
Puerto Escondido	inactiva o no reportada	INDRHI	meteoclimática	pública	-71.57233	18.32106
La Peñita	inactiva o no reportada	INDRHI	meteoclimática	pública	-70.27561	18.45551
Sabana de la Mar	regular	INDRHI	meteoclimática	pública	-69.38908	19.04605
Medina	activa o bueno	INDRHI	meteoclimática	pública	-70.14450	18.53523
Cítricos Dominicanos	inactiva o no reportada	INDRHI	meteoclimática	pública	-70.24028	18.73897
Engombe	regular	INDRHI	meteoclimática	pública	-70.00194	18.44941
Hato Mayor	inactiva o no reportada	INDRHI	meteoclimática	pública	-69.25001	18.75023
La Guama	inactiva o no reportada	INDRHI	meteoclimática	pública	-68.66119	18.80965
Palo de Caja	regular	INDRHI	meteoclimática	pública	-70.40006	18.53078
Presa Jiguez	regular	INDRHI	meteoclimática	pública	-70.38255	18.53328
Valdesia	regular	INDRHI	meteoclimática	pública	-70.28061	18.40856
Comate	inactiva o no reportada	INDRHI	meteoclimática	pública	-69.61664	18.74568
Pueblo Nuevo	inactiva o no reportada	INDRHI	meteoclimática	pública	-70.17477	18.43995
Triple Ozama	inactiva o no reportada	INDRHI	meteoclimática	pública	-69.94047	18.81030
Santana	inactiva o no reportada	INDRHI	meteoclimática	pública	-70.21687	18.27898
Los Jengibres	inactiva o no reportada	INDRHI	meteoclimática	pública	-70.04727	19.43913
Higuey	inactiva o no reportada	INDRHI	meteoclimática	pública	-68.69901	18.61992
El Memiso	activa o bueno	INDRHI	meteoclimática	pública	-70.57159	18.51421
Los Arroyos	activa o bueno	INDRHI	meteoclimática	pública	-70.52916	18.64196
Naranjal	inactiva o no reportada	INDRHI	meteoclimática	pública	-70.47712	18.54298
Trepada Alta	inactiva o no reportada	INDRHI	meteoclimática	pública	-69.45930	18.97636
El Seibo	inactiva o no reportada	INDRHI	meteoclimática	pública	-69.04444	18.76550
Jarabacoa	activa o bueno	INDRHI	meteoclimática	pública	-70.63896	19.13080
La Antona	activa o bueno	INDRHI	meteoclimática	pública	-71.40288	19.63358
Matagrande	activa o bueno	INDRHI	meteoclimática	pública	-70.98758	19.20857
Boca de Mao	inactiva o no reportada	INDRHI	meteoclimática	pública	-71.05147	19.58830
El Arroyaso	inactiva o no reportada	INDRHI	meteoclimática	pública	-70.53229	19.04496
El Naranjito	inactiva o no reportada	INDRHI	meteoclimática	pública	-71.34351	19.47442
Isa	inactiva o no reportada	INDRHI	meteoclimática	pública	-70.74591	19.44608
Quinigua	inactiva o no reportada	INDRHI	meteoclimática	pública	-70.77368	19.52663
Tavera	inactiva o no reportada	INDRHI	meteoclimática	pública	-70.71813	19.28357
Constanza	activa o bueno	INDRHI	meteoclimática	pública	-70.71674	18.91134
Vallejuelo	activa o bueno	INDRHI	meteoclimática	pública	-71.33949	18.65774
Km 11 Carretera Sánchez	inactiva o no reportada	INDRHI	meteoclimática	pública	-71.15053	18.75813
Peñón Barahona	inactiva o no reportada	INDRHI	meteoclimática	pública	-71.18722	18.29662
Angelina	activa o bueno	INDRHI	meteoclimática	pública	-70.22229	19.13080
Juma-Bonao	activa o bueno	INDRHI	meteoclimática	pública	-70.38618	18.90023
Barraquito	inactiva o no reportada	INDRHI	meteoclimática	pública	-69.78891	19.12899
El Aguacate	inactiva o no reportada	INDRHI	meteoclimática	pública	-69.75949	19.17220
Los Botados	inactiva o no reportada	INDRHI	meteoclimática	pública	-70.57674	18.87078
Blanco (Presa)	regular	INDRHI	meteoclimática	pública	-70.55812	18.88745
Cenovi en Santa Ana	regular	INDRHI	meteoclimática	pública	-70.34796	19.29467
El Limón	regular	INDRHI	meteoclimática	pública	-69.81945	19.15301
José Contreras	regular	INDRHI	meteoclimática	pública	-70.45006	19.46691
Presa Hatillo	regular	INDRHI	meteoclimática	pública	-70.25283	18.94745
Batey Excabacion	inactiva o no reportada	INDRHI	hidrométrica	pública	-69.40147	18.63061
Puerto Escondido	activa o bueno	INDRHI	hidrométrica	pública	-71.58066	18.32454
El Millo	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.34922	18.54715
Anamuya	inactiva o no reportada	INDRHI	hidrométrica	pública	-68.64482	18.69947
Rinconcito	activa o bueno	INDRHI	hidrométrica	pública	-71.76762	18.95801
Sonador	activa o bueno	INDRHI	hidrométrica	pública	-71.57687	18.70537
Puertecito	activa o bueno	INDRHI	hidrométrica	pública	-71.50987	18.80093
Cajuilito	activa o bueno	INDRHI	hidrométrica	pública	-71.59966	19.08263
Ranchitos	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.61729	18.93246
Pozo Hondo	regular	INDRHI	hidrométrica	pública	-71.48591	19.01005
Pedro Santana	regular	INDRHI	hidrométrica	pública	-71.69024	19.10257
El Estrecho	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.94147	19.80553
Imbert	regular	INDRHI	hidrométrica	pública	-70.82714	19.29927
El Recodo	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.33945	18.37361
Los Jengibres	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.04242	19.44386
La Pinta	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.55150	19.64413
Manabao	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.78760	19.07700
La Espensa	regular	INDRHI	hidrométrica	pública	-71.55428	19.57802
Santa Lucía	inactiva o no reportada	INDRHI	hidrométrica	pública	-68.93690	18.68206
El Mamey	activa o bueno	INDRHI	hidrométrica	pública	-68.75111	18.66878
Los Corozos	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.12344	18.52466
Don Miguel	activa o bueno	INDRHI	hidrométrica	pública	-71.68351	19.49365
Carbonera	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.69400	19.59052
La Aduana	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.71141	19.54847
Cinta Negra	inactiva o no reportada	INDRHI	hidrométrica	pública	-69.96343	19.28035
Villa Nizao	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.18576	18.01662
Bocaina	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.47316	18.67280
Jiguey	regular	INDRHI	hidrométrica	pública	-70.37745	18.54599
Aguacate	regular	INDRHI	hidrométrica	pública	-70.32149	18.49588
El Salto	activa o bueno	INDRHI	hidrométrica	pública	-69.60801	18.83988
El Higuero	inactiva o no reportada	INDRHI	hidrométrica	pública	-69.99933	18.61658
Don Juan	inactiva o no reportada	INDRHI	hidrométrica	pública	-69.95100	18.81812
Los Hidalgos	inactiva o no reportada	INDRHI	hidrométrica	pública	-69.66440	18.87713
El Cacique	regular	INDRHI	hidrométrica	pública	-69.86123	18.81733
Valle de Samaná	inactiva o no reportada	INDRHI	hidrométrica	pública	-69.33747	19.26323
San Rafael I	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.13787	18.02911
San Rafael II	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.14120	18.03241
San Rafael	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.14116	18.03042
Paso al Medio	inactiva o no reportada	INDRHI	hidrométrica	pública	-69.17193	18.56742
Santa Cruz	activa o bueno	INDRHI	hidrométrica	pública	-71.52177	19.62386
Puente Guayubin	activa o bueno	INDRHI	hidrométrica	pública	-71.39816	19.66163
Sabaneta	activa o bueno	INDRHI	hidrométrica	pública	-70.97864	19.54565
La Fortaleza	activa o bueno	INDRHI	hidrométrica	pública	-70.83344	19.04703
Puente San Rafael	activa o bueno	INDRHI	hidrométrica	pública	-71.05942	19.58971
Peña Ranchadero	activa o bueno	INDRHI	hidrométrica	pública	-71.32977	19.65942
Inoa	activa o bueno	INDRHI	hidrométrica	pública	-70.95897	19.32357
Boma	activa o bueno	INDRHI	hidrométrica	pública	-70.66224	19.16861
Jinamagao	activa o bueno	INDRHI	hidrométrica	pública	-70.98202	19.54277
Jarabacoa	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.67840	19.08357
Palo verde	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.55546	19.75778
Mata Grande	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.99611	19.19532
Cana Chapetón	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.25967	19.59789
Bejucal	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.73507	19.18246
Rincón	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.38663	19.52613
La Antona	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.40289	19.63025
Inaje	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.46927	19.42135
Jicomé	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.94842	19.62802
La Gorra	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.49482	19.51830
Bulla	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.07707	19.41982
Chorrera	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.08481	19.46413
Martínez	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.07953	19.47052
Mata de Jobo	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.33121	19.46497
Manabao	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.78987	19.06973
Paso Bajito	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.66399	19.01495
Bao	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.79702	19.30579
Sabana Iglesia	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.74619	19.31413
El Puente	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.76563	19.45218
Guanajuma	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.74535	19.28108
Jamamú	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.93344	19.19440
Jamamú	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.93344	19.19440
Casa de Máquina	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.59787	19.09892
Hato Viejo	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.63414	19.13474
Pinalito	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.75897	19.29524
Potrero	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.96008	19.46969
Paso de la Perra	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.83966	19.07752
Pinar Quemado	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.66513	19.09502
Los Velazquitos	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.67757	19.21274
Tavera	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.71535	19.26135
Baitoa	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.70896	19.32774
Baitoa	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.70896	19.32774
López	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.71118	19.33413
López	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.71118	19.33413
Las Charcas	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.71396	19.40857
Pte. Viejo Santiago	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.71174	19.43107
Pte. Viejo Santiago	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.71174	19.43107
Santiago	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.70313	19.44135
Santiago	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.70313	19.44135
Hato Yaque	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.79619	19.50580
Cañeo	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.98314	19.54247
Jicome	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.14731	19.62497
La Caña	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.52538	19.62247
Dique Juan Gómez	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.41705	19.68913
El Higuero	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.83769	19.24795
Redondo	regular	INDRHI	hidrométrica	pública	-70.58389	19.02863
Aguas Caliente	regular	INDRHI	hidrométrica	pública	-70.89814	19.24524
Pinalito	regular	INDRHI	hidrométrica	pública	-70.75897	19.29524
Inoa	regular	INDRHI	hidrométrica	pública	-70.98300	19.35222
Paso de la Palma	regular	INDRHI	hidrométrica	pública	-71.51604	19.57455
Los Pilones	regular	INDRHI	hidrométrica	pública	-70.77352	19.25579
El Cerrazo	regular	INDRHI	hidrométrica	pública	-70.72710	19.23592
Palomino	activa o bueno	INDRHI	hidrométrica	pública	-70.96944	18.80512
Paso de Lima	activa o bueno	INDRHI	hidrométrica	pública	-71.29826	19.02472
Vallejuelo	activa o bueno	INDRHI	hidrométrica	pública	-71.31721	18.65660
Hato Viejo	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.25538	18.79707
Palo Alto	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.17364	18.29928
Conuquito	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.15148	18.41161
El Rodeo	inactiva o no reportada	INDRHI	hidrométrica	pública	-71.26760	19.00551
La Coja	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.81868	18.76856
La Guama	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.81296	18.72572
Los Guiros	regular	INDRHI	hidrométrica	pública	-71.00095	18.51058
El Cacheo	regular	INDRHI	hidrométrica	pública	-71.11519	18.81385
Paso de Lima	regular	INDRHI	hidrométrica	pública	-71.29826	19.02472
Yásica Arriba	activa o bueno	INDRHI	hidrométrica	pública	-70.59327	19.63365
Palma Herrera	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.41089	19.55469
Sonador	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.53943	19.59870
Jamao	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.44828	19.63171
Los Brazos	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.43089	19.66052
Santa Ana	activa o bueno	INDRHI	hidrométrica	pública	-70.37064	19.25907
Yuna	activa o bueno	INDRHI	hidrométrica	pública	-70.25278	18.94722
Los Arroces	activa o bueno	INDRHI	hidrométrica	pública	-70.16638	18.97866
Blanco (Presa)	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.52139	18.88278
El Torito	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.45590	18.75634
Los Plátanos	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.31890	18.78632
Blanco	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.52146	18.88300
Rodeo	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.56840	18.88439
Tireo	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.63285	18.91912
Los Guázaros	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.52646	18.85579
Tireo	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.56840	18.88439
El Pino de Yuna	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.45618	18.79550
Piedra los Veganos	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.47173	18.81995
Hatillo	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.25283	18.94745
Arroyón Arriba	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.52827	18.85972
Arroyón Abajo	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.52354	18.86066
Jeremías	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.48597	19.22554
Jamo	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.48885	19.22868
Cañabón	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.41887	19.05395
Cevicos	inactiva o no reportada	INDRHI	hidrométrica	pública	-69.98432	19.00617
Sabana Grande	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.06621	19.03554
Jayaco	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.43936	19.00590
Blanco	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.36810	18.88717
Naranjal	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.51450	19.32830
Maimón	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.28490	18.89937
Los Guázaros	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.52257	18.85886
El Meche	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.53291	18.84794
Rincón de Yuboa	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.36157	18.78331
Los Quemados	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.47314	18.87916
Bonao	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.40762	18.95916
Los Corosos	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.20089	19.05496
La Jagua	inactiva o no reportada	INDRHI	hidrométrica	pública	-69.70796	19.14882
La Bija	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.12723	19.14948
Bayacanes	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.58153	19.23565
Cenoví, Los Limones	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.34906	19.31966
El pino de Yuna	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.29723	18.89472
Abadesa	inactiva o no reportada	INDRHI	hidrométrica	pública	-69.92643	19.01863
Los Plátanos	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.11228	19.12214
Ranchito	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.38872	19.17508
Los Tres Pasos	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.07879	18.97640
La Cabirma	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.11522	18.97470
Pinalito	inactiva o no reportada	INDRHI	hidrométrica	pública	-70.28923	18.89662
La Boca	regular	INDRHI	hidrométrica	pública	-70.30209	19.24578
Jaya	regular	INDRHI	hidrométrica	pública	-70.30250	19.24556
Rincón	regular	INDRHI	hidrométrica	pública	-70.40695	19.10509
Villa Riva	regular	INDRHI	hidrométrica	pública	-69.90444	19.17440
El Limón	regular	INDRHI	hidrométrica	pública	-69.81921	19.15301
BAIGUAQUEPILONES LOS	activa o bueno	INDRHI	hidrométrica	pública	-70.87945	19.34889
BAIGUAQUEPILONES LOS	activa o bueno	INDRHI	hidrométrica	pública	-70.87945	19.34889
Jardín Botánico de Santo Domingo	activa o bueno	INTEC	meteoclimática	privada	-69.95306	18.49500
San Pedro de Macorís	activa o bueno	INTEC	meteoclimática	privada	-69.28655	18.47433
Las Terrenas	activa o bueno	INTEC	meteoclimática	privada	-69.55257	19.03248
Jardín Botánico de Santiago	activa o bueno	INTEC	meteoclimática	privada	-70.69772	19.49751
Colegio Quisqueya, Santo Domingo	activa o bueno	INTEC	meteoclimática	privada	-69.95045	18.45637
Campus PUCCM	activa o bueno	INTEC	meteoclimática	privada	-70.68242	19.44092
ARPT. ARROYO BARRIL	activa o bueno	ONAMET	meteoclimática	pública	-69.43144	19.20050
ARPT. GREGORIO LUPERÓN	activa o bueno	ONAMET	meteoclimática	pública	-70.56316	19.75416
ARPT. MARÍA MONTEZ	activa o bueno	ONAMET	meteoclimática	pública	-71.12288	18.24861
BAYAGUANA	activa o bueno	ONAMET	meteoclimática	pública	-69.63083	18.74216
CABO ENGAÑO	activa o bueno	ONAMET	meteoclimática	pública	-68.32550	18.61678
CABRERA	activa o bueno	ONAMET	meteoclimática	pública	-69.90633	19.64436
CONSTANZA	activa o bueno	ONAMET	meteoclimática	pública	-70.71619	18.91173
HONDO VALLE	activa o bueno	ONAMET	meteoclimática	pública	-71.69364	18.71213
JARABACOA	activa o bueno	ONAMET	meteoclimática	pública	-70.60572	19.12551
JIMANÍ	activa o bueno	ONAMET	meteoclimática	pública	-71.85301	18.49278
ARPT. DE LA ROMANA	activa o bueno	ONAMET	meteoclimática	pública	-68.90925	18.44853
LA VEGA	activa o bueno	ONAMET	meteoclimática	pública	-70.53108	19.22275
LA VICTORIA	activa o bueno	ONAMET	meteoclimática	pública	-69.83956	18.59779
ARPT. DE LAS AMÉRICAS	activa o bueno	ONAMET	meteoclimática	pública	-69.67959	18.43306
LOS LLANOS	activa o bueno	ONAMET	meteoclimática	pública	-69.49096	18.62499
LOYOLA (SCR )	activa o bueno	ONAMET	meteoclimática	pública	-70.11273	18.41193
MOCA	activa o bueno	ONAMET	meteoclimática	pública	-70.53213	19.38811
MONCIÓN	activa o bueno	ONAMET	meteoclimática	pública	-71.16163	19.40988
MONTE CRISTI	activa o bueno	ONAMET	meteoclimática	pública	-71.65350	19.84988
OVIEDO	activa o bueno	ONAMET	meteoclimática	pública	-71.40403	17.80205
PADRE LAS CASAS	activa o bueno	ONAMET	meteoclimática	pública	-70.94115	18.73026
PERALTA	activa o bueno	ONAMET	meteoclimática	pública	-70.75596	18.57806
POLO	activa o bueno	ONAMET	meteoclimática	pública	-71.27231	18.08746
RANCHO ARRIBA	activa o bueno	ONAMET	meteoclimática	pública	-70.44973	18.69438
RÍO SAN JUAN	activa o bueno	ONAMET	meteoclimática	pública	-70.07475	19.63251
SABANA DE LA MAR	activa o bueno	ONAMET	meteoclimática	pública	-69.38883	19.05267
SALCEDO	activa o bueno	ONAMET	meteoclimática	pública	-70.41498	19.37531
SAMANÁ	activa o bueno	ONAMET	meteoclimática	pública	-69.33286	19.20438
SAN JOSÉ DE OCOA	activa o bueno	ONAMET	meteoclimática	pública	-70.50950	18.54200
SAN RAFAEL DEL YUMA	activa o bueno	ONAMET	meteoclimática	pública	-68.67416	18.42610
SÁNCHEZ	activa o bueno	ONAMET	meteoclimática	pública	-69.61311	19.23314
SANTIAGO	activa o bueno	ONAMET	meteoclimática	pública	-70.69806	19.50054
SANTO DOMINGO ESTE	activa o bueno	ONAMET	meteoclimática	pública	-69.87053	18.47339
VILLA RIVA	activa o bueno	ONAMET	meteoclimática	pública	-69.91483	19.17993
VILLA VÁSQUEZ	activa o bueno	ONAMET	meteoclimática	pública	-71.45730	19.74207
VILLA ALTAGRACIA	activa o bueno	ONAMET	meteoclimática	pública	-70.17060	18.67582
ALTAMIRA	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.83400	19.67516
ANGELINA	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.21977	19.13138
ARPT. DE PUNTA CANA	inactiva o no reportada	ONAMET	meteoclimática	pública	-68.35941	18.56400
ARPT. DEL CIBAO	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.59781	19.40305
ARPT. EL CATEY	inactiva o no reportada	ONAMET	meteoclimática	pública	-69.73366	19.26700
ARPT. LA ISABELA	inactiva o no reportada	ONAMET	meteoclimática	pública	-69.98158	18.57696
ARROYO LORO	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.27818	18.81403
BANÍ	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.31014	18.27867
BÁNICA	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.69836	19.08156
BASE AÉREA DE SAN ISIDRO	inactiva o no reportada	ONAMET	meteoclimática	pública	-69.76989	18.49944
BOHECHÍO	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.99150	18.77747
BONAO	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.40500	18.93317
CABRAL	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.21920	18.24784
COTUÍ	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.15273	19.04825
DAJABÓN	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.70436	19.54325
DUVERGÉ	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.51814	18.37666
EL CERCADO	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.51884	18.72556
EL SEIBO	inactiva o no reportada	ONAMET	meteoclimática	pública	-69.03995	18.77438
ELÍAS PIÑA	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.70133	18.87600
ENRIQUILLO	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.24173	17.89165
GASPAR HERNÁNDEZ	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.27311	19.62383
GURABO	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.66713	19.49794
HATO MAYOR	inactiva o no reportada	ONAMET	meteoclimática	pública	-69.25557	18.77005
HIGÜEY	inactiva o no reportada	ONAMET	meteoclimática	pública	-68.70098	18.62045
IMBERT	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.83042	19.75626
JARDÍN BOTÁNICO NACIONAL	inactiva o no reportada	ONAMET	meteoclimática	pública	-69.95293	18.49489
JUMA BONAO	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.38549	18.90047
LA VEGA - IATESA	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.55246	19.25230
LA VEGA - M. A.	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.49850	19.19920
LAS MATAS DE FARFÁN	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.52819	18.87478
LAS MATAS DE SANTA CRUZ	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.50686	19.65356
LIMÓN DEL YUNA	inactiva o no reportada	ONAMET	meteoclimática	pública	-69.78488	19.11985
LOMA DE CABRERA	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.61228	19.41925
LUPERÓN	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.96424	19.88717
MAO	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.04966	19.58833
MICHES	inactiva o no reportada	ONAMET	meteoclimática	pública	-69.04221	18.98099
MONTE PLATA	inactiva o no reportada	ONAMET	meteoclimática	pública	-69.78160	18.80243
NAVARRETE	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.86739	19.56062
NEIBA	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.41889	18.48347
PARQUE MIRADOR NORTE	inactiva o no reportada	ONAMET	meteoclimática	pública	-69.93112	18.52763
PARQUE MIRADOR SUR	inactiva o no reportada	ONAMET	meteoclimática	pública	-69.93416	18.45298
RADIOSONDA	inactiva o no reportada	ONAMET	meteoclimática	pública	-69.87016	18.47416
RESTAURACIÓN	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.68913	19.31497
SABANA GRANDE DE BOYÁ	inactiva o no reportada	ONAMET	meteoclimática	pública	-69.80099	18.94384
SAN CRISTÓBAL	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.09011	18.40710
SANTIAGO RODRÍGUEZ	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.33899	19.47712
TÁBARA ABAJO	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.86903	18.46438
TAMBORIL	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.61141	19.48342
VILLA ISABELA	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.05694	19.82111
VILLA LOS ALMÁCIGOS	inactiva o no reportada	ONAMET	meteoclimática	pública	-71.44328	19.41080
YAMASÁ	inactiva o no reportada	ONAMET	meteoclimática	pública	-70.02175	18.77271
Río Ozama en barrera de salinidad	activa o bueno	SGN	hidrométrica	pública	-69.78653	18.57151
Río Haina en Manoguayabo	activa o bueno	SGN	hidrométrica	pública	-70.01934	18.47068
Estación Japonesa	activa o bueno	Universidad ISA	meteoclimática	privada	-70.74828	19.44320